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'^ ' Abstract 

^ , The main part of this review is devoted to the comprehensive description of coher- 

ent radiation by nuclear spins. The theory of nuclear spin superradiance is developed 
and the experimental observations of this phenomenon are considered. The intriguing 
^ . problem of how coherence develops from initially incoherent quantum fluctuations is 

Cn \ analysed. All main types of coherent radiation by nuclear spins are discussed, which 

are: free nuclear induction, collective induction, maser generation, pure superradi- 
\C) • ance, triggered superradiance, pulsing superradiance, punctuated superradiance, and 

Cp I induced emission. The influence of electron-nuclear hyperfine interactions and the role 

f— ^ ' of magnetic anisotropy are studied. Conditions for realizing spin superradiance by mag- 

netic molecules are investigated. The possibility of nuclear matter lasing, accompanied 
by pion or dibaryon radiation, is briefly touched. 



1 Introduction 

Nuclei can radiate in different ways. For example, they emit gamma radiation in the process 
of changing internal quantum states. An ensemble of nuclei, emitting such hard photons, 
could form a coherent source, called gamma-ray laser or just gamma laser. However, the 
problem of gamma lasers remains a challenging but frustrating field of research, with not a 
great progress in theory and practically no experimental achievements [1-3]. 

Contrary to this, there exists a type of coherent nuclear radiation that is well documented 
experimentally and for which a detailed microscopic theory has recently been developed. 
This is nuclear spin radiation. The main part of the present review is just devoted to the 
phenomenon of coherent radiation by nuclear spins. We, first, explain in simple words the 
physics of this effect and survey the basic experiments where it was observed. Then we pass 
to developing a comprehensive microscopic theory of strong nonlinear dynamics of nuclear 
spins. A special attention is paid to a very intriguing problem of how coherence arises from an 
incoherent quantum motion of randomly fluctuating spins. We describe all principal regimes 
of spin superradiance and study the influence of the hyperfine electron-nuclear coupling. 
We analyse the role of the single-ion magnetic anisotropy that is so important in the spin 
radiation of magnetic molecules. Finally, we note that not only photons can be emitted 
coherently by nuclei. Thus, excited nuclear matter can produce coherent radiation of gluons, 
pions, and dibaryons. Different types of coherent nuclear radiation can find a variety of 
applications, some of which are discussed in the review. The content of the latter is as 
follows: 

1. Introduction 

2. Nuclear Spin Superradiance 

3. Ensemble of Nuclear Spins 

4. Nuclear Spin Waves 

5. Scale Separation Approach 

6. Incoherent Quantum Stage 

6.1. Radiation by Magnetic Dipoles 

6.2. Resonator Nyquist Noise 

6.3. Local Spin Fluctuations 

7. Regimes of Coherent Radiation 

7.1. Transient Spin Superradiance 

7.2. Pulsing Spin Superradiance 

7.3. Induced Coherent Emission 

8. Electron-Nuclear Hyperfine Coupling 

9. Enhanced Nuclear Radiation 



10. Superradiance by Magnetic Molecules 

11. Pion and Dibaryon Radiation 

12. Conclusion 

2 Nuclear Spin Superradiance 

Atomic systems, radiating at optical frequencies, exhibit a number of coherent effects [4,5]. 
The majority of the latter have their counterparts in spin systems, generating radiation at 
radio-frequencies. One of the most interesting coherent effects is superradiance that may 
occur in both atomic and spin systems. The possibility of superradiance in atomic systems 
was predicted by Dicke [6] and the feasibility of organizing coherent motion of spins was 
discussed by Bloembergen and Pound [7]. The modern status of nuclear spin superradiance 
is presented in review [8]. 

By definition, superradiance is the process of coherent spontaneous radiation. Being 
spontaneous, it is assumed to be self-organized. In general, a coherent motion of spins can 
be realized by means of a strong external magnetic field. But this would result in nuclear 
induction, which, though is a coherent process, however is not superradiance. To realize the 
latter, one, first, needs to prepare the spin ensemble in a nonequilibrium state. For this, it is 
possible to polarize spins in one direction and after that to place them in an external magnetic 
field of opposite direction. Such a system of spins would be analogous to an ensemble of 
inverted atoms. Several other similarities between atomic and spin systems are discussed in 
reviews [8,9]. Two important problems related to spin assemblies are: What is the cause 
provoking spins to start their initial motion, after they have been inverted? And what is 
the mechanism collectivizing the following spin motion, making the latter coherent? The 
answer to the second question was proposed by Bloembergen and Pound [7] who suggested 
to place the spin sample into an electric coil being a part of an electric circuit with a 
natural frequency tuned close to the Zeeman transition frequency of spins. This coupling 
with a resonant electric circuit would produce a feedback field collectivizing the spin motion. 
The answer to the first question, what is the cause starting the spin motion, has not been 
understood for many years. This puzzle was solved recently [10] and will be considered in 
detail in the following sections. In brief, there are quantum fluctuations of spins that trigger 
the initial spin motion. 

One more question would be how to measure the radiation generated by moving spins? 
The magnetodipole radiation, resulting from this motion, even being completely coherent, 
would have the intensity of order 

m - ^ (^0 ir co^ N' , 

where c is the light velocity; /iq = ^7n, with 7„ being the nucleus gyromagnetic ratio; / is the 
nuclear spin; ujq is the Zeeman transition frequency; A^ is the total number of spins taking 
part in the coherent radiation. Accepting the typical values /xq ~ 10^^^ erg/G, ujq ~ 10* 
Hz, and A^ ~ 10^^, we would have I{t) ~ 10~^ W, where 1 W= 10'^ erg/s is a watt. This 
is a rather small quantity of the intensity /(t), which would be different to notice. But this 



radiation induces electric current in the coil, with the power of current P(t) being essentially 
larger than the intensity I{t). Roughly speaking, 
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i{t) K ' 

where A is the radiation wavelength and Vc is the coil volume. This ratio, for A ~ 10^ — 10^ 
cm and Vc ~ 10 cm^, is of order 10^ — 10^, which makes P{t) an easily measurable quantity. 
Thus, the magnetodipole radiation of a spin sample is quite weak and practically does not 
propagate into free space, but mainly is taken up by a resonant coil surrounding the sample 
[11]. 

One should not confuse superradiance with other types of coherent radiation. By defi- 
nition, superradiance is spontaneous collective emission. The term spontaneous implies that 
the process is self- organized but not induced by external fields. And the word collective 
means that the radiation characteristics are essentially influenced by collective effects. For 
instance, the radiation intensity of A^ radiators is proportional to A^", with a > 1, while the 
radiated pulse is short, with the radiation time proportional to A^^"". Depending on the 
relation between characteristic times, there can arise different kinds of coherent radiation. 
Among these typical times, the most important are: The crossover time tc, separating the 
quantum incoherent stage of spin motion from their coherent rotation; the pulse time Tp 
of an emitted pulse; and the dephasing time T2. In the case, when there are no external 
transverse fields, one can distinguish seven coherent radiation regimes: 

(1) Free induction {tc = 0, Tp> T2). 

Coherence of radiation here is not self-organized but imposed upon the system by an 
external field. This is the standard nuclear free decay induced by an initial coherent pulse 
or, equivalently, by an essential transverse polarization [12]. 

(2) Collective induction (tc = 0, Tp < T2). 

The shortening of the radiation damping is due to the coupling with a resonator. However, 
the process is mainly governed not by a self-organization but by a strong external coherent 
field or by an initial transverse polarization [7]. 

(3) Maser generation {tc > 0, Tp > T2). 

There are no external coherent fields, though some incoherent nonresonant pumping may 
exist. Self-organization becomes crucial. But the emitted pulses are not sufficiently narrow, 
implying that self-organization is not yet high. Although this is a spontaneous coherent 
process, it is not yet a genuine superradiance. Such a type of coherent maser generation 
was observed in many experiments [13-16]. The physics of this process is analogous to the 
known laser generation [17-19]. 

(4) Pure superradiance {tc > 0, Tp < T2). 

This is a purely self-organized process, developing from an incoherent chaotic stage to 
a highly correlated spin motion. Though some of the features of spin superradiance are 
similar to those of atomic superradiance [4,20-22], there are also several principal differences. 
Experiments on and theory of spin superradiance will be described in detail below. 

(5) Triggered superradiance {tc > 0, Tp < T2). 



In this process, self-organization is crucially important, but there is as well a weak ex- 
ternal resonant field triggering and infiuencing the behaviour of spins. There are similarities 
with triggered superradiance in optics [20]. 

(6) Pulsing superradiance (tc > 0, Tp < T2). 

This regime differs from pure superradiance by the existence of a series of superradiant 
pulses, instead of a single superradiant burst. To realize such a pulsing superradiance, it is 
necessary to maintain a sufficiently high level of spin polarization by means of an incoherent 
pumping. The latter can be done, e.g., by dynamic nuclear polarization [12,23,24]. 

(7) Punctuated superradiance (tc > 0, Tp < T2). 

Here, similar to pulsing superradiance, there is a series of superradiant pulses; however 
this is achieved not by supporting nuclear polarization but by means of resonant external 
fields or forces. All regimes of spin superradiance will be thoroughly considered in the 
following sections. 

From the above classification it follows that superradiance, in addition to being a coherent 
and spontaneous emission, should correspond to short pulses with a finite crossover time, so 
that 

tc> , Tp <T2 . 

Thus, the correct and full definition of superradiance would be: 

Superradiance is a coherent spontaneous emission by an ensemble of radiators of short 
pulses peaked at a time that is larger than the crossover time, with the duration of each pulse 
shorter than the dephasing time. 

This definition suits for any kind of superradiance, whether it is atomic or spin super- 
radiance. In what follows, we concentrate on spin superradiance produced by an ensemble 
of nuclear spins. Initially, we shall give a brief survey of experiments observing nuclear spin 
superradiance. 

Pure spin superradiance by nuclei was first observed in Dubna [25,26] and confirmed in 
St. Petersburg [27,28]. These experiments were accomplished with propanediol C3H8O2. 
This is a material rich of protons, with the density p/^ ^ 4 x 10^^ cm~^. The proton spins 
were polarized, by dynamic nuclear polarization, in an external magnetic field i?o ~ 1 T, 
which corresponds to the Zeeman frequency uj^ ~ 10^ Hz. The sample was refrigerated 
to low temperatures T ~ 0.1 K, at which the nuclear spin-lattice relaxation was strongly 
suppressed, with the longitudinal relaxation time Ti ~ 10^ s. The transverse dephasing time, 
due to dipole spin-spin interactions, was T2 ~ 10^^ s. The coupled resonant electric circuit 
had a quality factor Q ~ 100 and a ringing time r ~ 10^^ s. 

Similar experiments were accomplished in Bonn [29] with butanol C4H9OH and ammonia 
NH3. In these materials, the proton density pn ~ 10^^ cm^^. The characteristic Zeeman 
frequency was uq ~ 10^ Hz. The materials were kept at low temperature, when the spin- 
lattice relaxation is suppressed, with the longitudinal relaxation time Ti ~ 10^ s. The 
transverse dephasing time was T2 ~ 10~^ s. The quality factor of the resonator electric 
circuit was Q ~ 30, and the ringing time r ~ 5 x 10~^ s. 

Such proton-rich materials are widely used as targets in studying the scattering of particle 
beams from accelerators. It would, of course, be possible to involve for experiments on spin 



superradiance other materials employed as targets in high-energy scattering studies. For 
instance, a good candidate would be pentanol C5H12O, whose proton density is p// ~ 6 x 10^^ 
cm ~^. 

Pulsing spin superradiance was explored in Ziirich [30-33] on the basis of nuclei of ^^Al 
inside the ruby crystal AI2O3, where the density pAi ~ 4 x 10^^ cm^^. The nuclei ^'''Al 
have spins / = 5/2. The crystal was oriented in an external magnetic field i?o ~ 1 T 
so that the fully resolved structure of the five Am = ±1 transition lines could be clearly 
seen. When a resonant circuit is tuned to a selected transition line, then ^^Al spins form 
an effective two-level system. In experiments [30-33], the circuit was usually tuned to the 
central { — 1,1} line, with a transition frequency uq ~ 10® Hz. At low temperatures around 
T ~ 1 K, the spin-lattice relaxation time was Ti ~ 10^ s and the transverse dephasing time 
was T2 ~ 10"^ s. The resonant circuit had the quality factor Q ~ 100 and ringing time 
r ~ 10^^ s. The inversion of spin polarization was permanently supported by means of 
dynamic nuclear polarization with the pumping rate FJ between 0.01 and 10 s~^. 

A typical experimental setup employed for observing spin superradiance [25,26] is shown 
in Fig. 1. The detected superradiance pulse has the form illustrated in Fig. 2 for two 
different initial spin polarizations. The higher is the prepared inversion, the stronger is the 
superradiant burst. 

The influence of the passive resonant circuit, coupled to a nuclear spin system, was also 
studied in nuclear spin echo experiments [34-36]. In the ferrite Lio.5Fe2.5O4, the spins of 
nuclei ^^Fe were tuned to a resonant circuit with a frequency about 7 x 10^ Hz. In the 
compound Co2MnSi, the working substance was the nuclei ^^Co, whose spins were coupled 
to a circuit of frequency 1.4 x 10® Hz. The experiments [34-36] showed that the presence of 
the resonant circuit can enhance the intensity of the nuclear spin echo signal by up to 40%, 
as compared to the case without coupling to such a resonant circuit [37]. 

The theoretical description of spin dynamics in the presence of a resonant circuit has been 
commonly done on the basis of the Bloch plus Kirchhoff equations [27-33,38-40]. Analytical 
solutions were usually obtained by involving the adiabatic approximation. However, the 
latter approximation is only valid for dynamics close to stationary and it does not correctly 
describe transient effects, such as superradiant bursts. What is more, the Bloch equations 
presuppose the existence of coherent motion from the beginning, since these are classical 
equations. Such equations principally are not able to depict self-organized coherent phe- 
nomena, as pure superradiance, which develops from an initially incoherent quantum stage. 
This principal defect of Bloch equations cannot be overcome even if one employs more elab- 
orate approximations [41,42] based on the averaging technique [43]. As was first shown in 
Ref. [10], for the correct description of pure spin superradiance it is crucially important to 
take account of stochastic local spin fluctuations. This conclusion was also confirmed by 
numerical calculations [44-46]. 

Computer simulations of spin superradiance has been accomplished [47-51], investigating 
different regimes of coherent spin motion. A typical result of simulations is shown in Fig. 3. 
The weak points of such simulations are: First, one has to deal with a rather limited number 
of spins, say 10^ or 10^, because of which the obtained solutions give only a qualitative 
picture. Second, in such simulations, spins are treated as classical vectors, hence, quantum 
fluctuations, so important at the initial stage of motion, are not properly considered. Third, 
numerical results do not always present clear physical explanations of the studied processes. 
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This is especially so, when there are tens of parameters to be varied, while in reality only 
some combinations of these parameters are important. 

To give a thorough description of nuclear spin superradiance, it is necessary to develop a 
good microscopic theory of spin dynamics in a sample coupled to a resonant electric circuit. 
The feedback field created by the resonator is crucial for organizing the coherent motion 
of spins, moving from a strongly nonequilibrium state. Note that for spin motion close 
to equilibrium, as in nuclear magnetic resonance experiments, the feedback field is not of 
much significance [52]. For describing such a purely self-organized process as pure spin 
superradiance, it is vital to take into account local quantum spin fluctuations, which trigger 
the initial spin motion and define the delay time of the superradiant burst [10]. 

A comprehensive microscopic theory of spin superradiance, based on realistic spin Hamil- 
tonians [12,23,53,54] has been developed [10,55-60]. This theory for the first time has made 
it possible to discover the actual origin of pure spin superradiance [10,55-57], to give an 
accurate description of various regimes of nonlinear spin dynamics [55-59], to present an 
explicit picture of how coherence develops from chaotic spin fluctuations [10,60], to consider 
spin superradiance in different materials, such as polarized nuclear paramagnets [10,56-59], 
ferromagnets, ferrimagnets [61-63], and molecular magnets [60], to analyse the superradiant 
operation of spin masers [60,64,65], and to advance the feasibility of punctuated spin super- 
radiance [66]. Brief survey of this theory was done in reviews [8,9]. In the following sections, 
we shall give a detailed account of the theory and of several its applications. 

3 Ensemble of Nuclear Spins 

Consider a solid sample containing A^ nuclear spins /» enumerated by the index i = 1,2, . . . ,N. 
The Hamiltonian of the nuclear system can be written as 

h = T.^^ + It.h^,, (1) 

i i^j 

with Hi being related to individual spins and Hij corresponding to spin interactions. The 
single-spin term 

H^ = -fioB ■ h (2) 

is the Zeeman energy, where /xq = gifJ'W = film with gj being the Lande factor for a nucleus 
of spin J; fj,iy, nuclear magneton; 7„, nuclear gyromagnetic ratio. Generally, for nuclei, /Xq 
can be positive as well as negative. The total magnetic field 

B = Boe, + {B, + H)e^ (3) 

contains external longitudinal, Bq, and transverse, Bi, magnetic fields, and also a feedback 
field H of the resonant electric circuit. The pair terms 

H^, = E <^r ^f (4) 

a/3 

in the Hamiltonian (1) correspond to dipolar spin interactions through the dipolar tensor 

D?f = 4 {S^p - 3<, nj) , (5) 



V 



in which a, P = x,y,z and 



fij — Yijl , rijj _ 




The dipolar tensor enjoys the properties 









^i- 



3, (6) 

of which the second is exact and the first one is asymptotically exact for a macroscopic (in 
all directions) sample with a large number of spins A^ ^ 1. 

The single-spin term (2), with the total magnetic field (3), can be written as 

H, = -/ioi^o/r - ^ /io (5i + H) (/+ + ir) , (7) 

where Ij = H ± ilj are the ladder spin operators. Employing the latter, the interaction 
term (4) can be presented in the form 

H.J = a., {i: I! - \ltlj)+h,ltlt + Vl,l:rlT + 2c.,lt^ + 2 4^l7l^, (8) 

in which the notation 

«.. = ^^ , c^,^\ {D!f - ^Dr) , 6,^ ^ i (d- -D^-2^ D^f) (9) 

is introduced. From the first of Eqs. (6), it follows that 

E «ii = E biJ = J2 Cij = 0. (10) 

Jiy^i) Jiy^i) 3(.+i) 

Note also that the quantities (9) are symmetric, so that ajj = Ojj, 6^- = 6jj, and Cjj = c^j. 

In the total magnetic field (3), the external longitudinal and transverse fields Bq and i?i, 
respectively, are supposed to be prescribed. The resonator feedback field H is created by 
the electric current of the coil surrounding the sample. The orientation of the coordinate 
axes with respect to the latter is illustrated in Fig. 4. The electric circuit is characterized 
by resistance i?, inductance L, and capacity C. The spin sample is inserted into a coil of n 
turns, length /, and cross-section area Ac- The electric current in the circuit is described by 
the Kirchhoff equation 



4+«^+^/„'^(*')*'-^/-f- ("; 



in which Ej is an electromotive force and 

$ = — n AcT] rria: (12) 

c 

is a magnetic flux formed by the a;-component of the magnetization density 

^^ = y E < ^r > , (13) 



where rj ^ V/Vc is a filling factor, V is the sample volume, and the brackets < . . . > imply 
statistical averaging. 

The electric current, circulating over the coil, creates a magnetic field 

H^fl,. ,U) 

Let us employ the standard notation for the circuit natural frequency 

and the circuit damping 

where Q is the resonator quality factor and r is called the circuit ringing time. Define the 
reduced electromotive force 

'-^. (17) 

Then the Kirchhoff equation (11) can be transformed to the equation 

^ + 2^H + uj^ f'H{t')dt' = ^h-47rr]^ (18) 

at Jo at 

for the feedback magnetic field created by the coil. 

The feedback equation (18) can be presented in another equivalent form that proved to 
be very convenient for solving the evolution equations [55-60] . For this purpose, by involving 
the method of Laplace transforms and employing the transfer function 
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G(t) = [cos u't sinu't) e-^' [u;' = Ju^ - -f^) , (19) 



we find the integral representation 

H = I Git- t') hhit') - A-Kr^m^it')] dt' , (20) 

Jo 

where the overdot means, as usual, time derivative. 

4 Nuclear Spin Waves 

In order to understand the nature of the mechanisms triggering the spin motion, it is neces- 
sary to pay attention to the possibility of arising nuclear spin waves. The latter are known to 
exist in ferromagnets, where ferromagnetism is caused by electron spins [54]. To produce spin 
waves, nuclear spins should be somehow ordered, at least locally. The dipolar interactions 
(4), by themselves, usually are not able to order nuclear spins. However, in the presence of 



a sufficiently strong external magnetic field, typical of nuclear magnetic resonance experi- 
ments, some kind of the intermediate range ordering may appear among nuclear spins. This 
mid-range ordering can be observed by different NMR techniques [67]. 

The situation with nuclear spins, forming a paramagnetic matter, seems to be analogous 
to what happens in paramagnetic systems of electron spins. Therefore, it is useful to make 
a retrospective journey into the problem of spin waves in electron paramagnetic assemblies. 
This may better explain the physics of the similar effects in nuclear spin samples. 

The theory of electron spin waves in nonmagnetic metals was developed by Silin [68,69] 
on the basis of the Landau Fermi-liquid phenomenological picture [70,71]. Good surveys of 
the Silin theory can be found in [72-74]. Spin waves in nonmagnetic metals can exist only in 
the presence of an external magnetic field, and do not exist when this field is absent. To be 
well defined, spin waves require that the external magnetic field be sufficiently strong, so that 
the spin resonance frequency be larger than the spin-wave attenuation. Spin waves, in the 
field of about 10^ G, were observed in several paramagnetic metals [75-78], such as Na, K, 
Rb, Cs, and Al. Calculation of electro-magnetic fields transmitted through metallic films by 
spin waves involves the usage of boundary conditions defined by the type of spin scattering 
on the surface of metals. Several kinds of these conditions were employed in the theory of 
electron spin resonance [79-87]. The Dyson boundary condition [79] was used for studying 
the influence of the surface magnetic anisotropy on the amplitudes of excited spin waves and 
on the transition coefficient under spin- wave resonance [89-91]. Actually, experiments on 
spin resonance cannot distinguish between different boundary conditions [85]. Therefore, the 
usage of the simple Dyson condition for calculating the signals transmitted by spin waves 
[89-91] was justified. The surface impedance of a semi-bounded metal was expressed through 
the characteristics of infinite matter [90] . The possibility of exciting spin waves in a semi- 
bounded paramagnetic sample was advanced [92,93]. Thus, electron spin waves do exist in 
paramagnetic metals, provided that a sufficiently strong external magnetic field is applied. 
All this sets us thinking that an external magnetic field could also ensure the appearance of 
nuclear spin waves in paramagnetic nuclei assemblies. 

Let us consider a paramagnetic system of nuclear spins, described by the Hamiltonian 
(1). The spin operators satisfy the Heisenberg equations of motion 



ih 



dll 

dt 



It.H 



and the commutation relations 



r r 

« ' 3 



26ij I- , 






±% If . 



To write down the evolution equations in an explicit way, it is convenient to introduce the 
notation 



e = r E (2q, /| - i a,, ir + 26,, I] 



h 



Jif^i) 



(21) 
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for the local fields acting on a spin from the side of other spins. For short, we do not label 
the local fields (21) by an index i. Also, define the effective force 

/^-^/io(i?i + //)+e. (22) 

Then the evolution equations for the spin operators can be presented in the form 

dlr 



/+ /- + /+/). (23) 



dt 

dll^_l 
dt 2 

Denoting the statistical average of an operator I" as < If >, we may write the identity 

/f ^ < If > +5 /r , 51^^ It- <ir>: (24) 

where 61°' describes an operator deviation from the average value < /" >. Also, without 
the loss of generality, we may assume that the averages 

<I°> = <If> (25) 

do not depend on the indices i,j, thus, transferring this dependence to the deviation term 
61°. Then, because of Eqs. (10), the local fields (21) become 

eo = ^ E («^. ^^1 + 4 ^^i + c^.- ^^j 

^ = ^ E (2q, 61] - \ a,, 6If + 26,, 6^\ . (26) 

This better clarifies the meaning of the local fields .^o = <^Co) and ^ = 6^^ actually, corre- 
sponding to local field fluctuations. 

To show that an external magnetic field Bq can support the appearance of nuclear spin 
waves, let us consider a stationary situation, when Bq = const and Bi = H = 0, so that 
there is a nonzero polarization < If >^ 0, but < I^ >= 0. In such a case, I^ = 6lf. 
Equating in Eqs. (23) the terms linear with respect to spin deviations, we have 

|/r = ^.oSo/-+ </,->«, |^/.' = 0. (27) 

Since (5/f = const, we may set, by accepting the zero initial condition, that 6 If = 0. For the 
raising and lowering spin operators, we employ the Fourier transforms 

^f = E ^k exp {T^k ■ r,) , 4^ = ^ E ^f ^xp (±^k ■ r,) . (28) 

k j 

Strictly speaking, the latter are exactly valid only if all spins were located in sites of an 
ideal crystalline lattice. When the studied spin sample does not form such an ideal lattice, 
transformation (28) can be treated as approximate. In that case, the summation over spins 
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should be limited by a finite number of them inside a volume, with an effective size -^e//; 
where the matter can be regarded as approximately arranged in a lattice. The length Lg// 
is an effective linear size of uniformity. 
Introduce the transforms 

ak= Yl ^ij 6^P (~^'^ ■ ^ij) ' bk= Yl ^ij ^^P (~^^ ■ "^ij) ■ (29) 

Since aij = aji and it is real, a^ is also real. And, because of Eqs. (10), ao = Bq = 0. Then, 
using Eqs. (28) and (29), we reduce the first of Eqs. (27) to the equation 

— 4" = -ifik Ik + i^k Ik , (30) 

in which 

^. = I < /; > - ^ /io 5o , A, = ^ < Jf > . (31) 

The solution to Eq. (30) can be presented as 

4" = Uk e-*"'=* + vl e'"'=* , (32) 

with the spectrum of excitations 



^k = \IijI - |AfcP . (33) 

The uniform excitation, corresponding to A; = 0, yields 

^0 = T lAto Bq\ . (34) 

n 

The procedure of defining the excitation spectrum, as is used above, is analogous to that 
employed for ferromagnets [94]. 

The spectrum (33) can describe nuclear spin waves if it is positive. The requirement 
Uk > leads to the stability condition 

\2fioBo- < If > afcl > 4 |< /; > h\ . (35) 

From here, it is evident that in a paramagnet without a polarization, when < Jf >= 0, and 
without an external field, i.e. Bq = 0, there are no nuclear spin waves. But if there exists a 
polarization and an external field, that is, < /f >7^ and Bq ^ 0, then condition (35) can 
be valid for a sufficiently strong external field Bq. Therefore, nuclear spin waves are possible 
in a paramagnet placed in a magnetic field. 

In the long-wave limit, the spectrum (33) is given by the equation 

ul ^ u;l+ < /; > ^ E «^. (k ■ r.,)^ , (36) 

where /c — > 0. In the interval 

L;}f <^k<^a-\ (37) 
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where Le// is an effective length of uniformity and a is the mean distance between spins, 
one can hmit itself by the summation over nearest neighbours only. As an illustration, let us 
assume cubic symmetry, when there are six nearest neighbours. Then the tensors, defined 
in Eq. (9), are 

{a,,} = ^ {1, 1, 1, 1, -2, -2} , {h,} = ^ {-1, -1, 1, 1, 0, 0} , (38) 

while Cij = 0. The Fourier transforms (29) become 

ttk = — 5^ (cos k^a + cos kya — 2 cos kzO) , bk = ^ (cos k^a — cos kyu) , 

which, in the case of Eq. (37), gives 

a,^-^ {kl + kl-2kt) , 6,^M {kl-k^;) . (39) 



a ^ " ^ Aa 

In this manner, Eq. (36) results in the spectrum equation 

ul ^ ul+ <F^>l^(kl + kl- 2kl) . (40) 

The spectrum Uk, to be positive at the maximal wave vector k = a'~^, requires that 



ha^Uo 



< 1 . 



In the standard situation, when ^\/%a^ ~ 10^ s^^ and ujq ~ 10^ Hz, this inequality holds 
true. The attenuation of spin waves is of order T2 ~ fi^/ha^ ~ 10^ s~^, hence it is much less 
than uo ~ 10^ Hz. Thus, nuclear spin waves are well defined. 

It is necessary to stress the importance of an external magnetic field for the stability of 
nuclear spin waves. Really, setting i?o = in condition (35), we have \ak\ > 4|6fc|. For a 
cubic structure with the wave vectors in the interval (37), this yields 



k^ + ky zkt. 



>3 



k^ ky 



This inequality should be valid for all k = {kx, ky, kz} from the given interval. However, as 
is evident, the inequality does not hold if either kx = k = |k| or ky = k. Hence, these nuclear 
spin waves are not stable without an external magnetic field. 

One could remember that at very low temperatures, of order T ~ 10^^ K, nuclear 
spins interacting through dipolar forces can become magnetically ordered [23]. Then in a 
spin system with long-range magnetic order, there should appear spin waves, even without 
external fields. This, really, may happen, but only for those structures that are able to 
support nuclear spin waves. Generally, a stable structure is characterized by two types of 
stability, thermodynamic and dynamic [95]. The structure is thermodynamically stable, 
when its thermodynamic potential is extremal. For instance, when free energy is minimal. 
And the structure is dynamically stable if its spectrum of collective excitations is non- 
negative. This is also true for magnetic structures formed by dipolar spin systems. Only 
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those magnetically ordered structures will survive, which possess the minimal free energy 
and are stable against arising spin waves. 

As has been demonstrated above, nuclear spin waves may exist even when there is no 
long-range magnetic order caused by internal forces, but provided that there is a sufficiently 
strong external magnetic field. The latter can stabilize nuclear spin waves for any kind of 
crystalline symmetry. Moreover, the spin sample can be polycrystalline, being composed 
of many small crystals, or even it can be amorphous. Then not all space of the sample 
will be filled by a coherent system of spin waves. But the whole sample will be separated 
into regions, inside each of which spin waves are coherent, though are not coherent between 
different regions. The effective linear size of such regions, where spins can form coherent spin 
waves, Leff, can be called the uniformity length. If spin waves arise inside the spatial regions 
of volume L^rr which is much smaller than the total sample volume, then such waves may 
be termed local spin waves. The coherence length L^ff must be much larger than the mean 
interspin distance in order that spin waves could arise, though this length may be smaller 
than the linear size of the total sample. 

5 Scale Separation Approach 

To study nonstationary spin dynamics, we have to deal with the general equations of motion 
(23). In order to analyse these equations, we invoke the scale separation approach [9,42,56,57] 
consisting of three main steps: Randomization of local fields; Classification of relative quasi- 
invariants; and Method of stochastic averaging. 

Let us recall that there exist the operator constructions ^o and ^, playing the role of local 
fields (21) or, according to Eq. (26), of local field fluctuations. We may consider these local 
fields as a separate kind of operators. Since these fields describe local fiuctuations, they 
can be modelled by random variables [12,53]. In this way, we have two types of variables 
characterizing the system, spin operators I^ , I^ , If and random fields ^q, ^, C,*. The 
former are responsible for long-range global phenomena while the latter, for short-range 
local fiuctuations. Having defined two types of variables, we may introduce two sorts of 
averaging with respect to the corresponding variables. Thus, the statistical averaging over 
spin operators will be denoted by the single angle brackets < . . . > and the averaging over 
random local fields will be denoted by means of the double angle brackets <^ ... ^. This 
modelling of local fluctuating flelds by random variables is the randomization of local fields. 
To make the method self-consistent, it is necessary to deflne the stochastic averages over the 
random flelds. Treating the latter as white noise, we deflne 

< Ut) > = < m > = , < Ut) Ut') > = 2r3 s{t - 1') , 
« Ut) iit') » = « m iit') » = , « Cit) m » = 2r3 S{t - 1') , (41) 

where T^ is the width of dynamic broadening., which is the inhomogeneous broadening caused 
by the local fleld fluctuations. This broadening is similar to that arising in optical resonant 
systems because of dipole-dipole interactions through the exchange of transversely polarized 
photons [96]. 

Averaging over spin operators, we could employ the decoupling 

</rjf > = </r></f > {i^j). (42) 
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This looks as the mean-field approximation. However, we should remember that the aver- 
aging, denoted by the single brackets <...>, by definition does not involve the stochastic 
variables. Therefore the quantum correlations are not lost in the decoupling (42) but are 
preserved because of the dependence of the averages < /" > on the random variables .^o and 
^. This decoupling (42) may be termed the stochastic mean-field approximation [9]. More- 
over, after having agreed to treat the spin operators and the local fields (21) as different 
variables, we do not often need the explicit usage of the decoupling (42). This will be used 
only once in the definition of coherence intensity. 

Let us average the operator equations (23) over the spin operators, without touching 
the local random fields (21). And let us introduce the following definitions. The transition 
function 

1 ^ 
« = — - V < /r > (43) 

describes the average rotation of transverse spin components. The coherence intensity 

1 ^ 

shows the degree of coherence in the spin motion. And the spin polarization 

1 N 

defines the average polarization in the system. Note that, under the validity of the decoupling 
(42), and for A^ ^ 1, we have w = |-up, with u given in Eq. (43). 

Assume that the longitudinal magnetic field Bq is constant and directed so that 

/io 5o < . (46) 

For nuclei with /iq > 0, this means that Bq < 0. Conversely, if /io < 0, then Bq > 0. To 
allow for the infiuence of lattice, account must be taken of spin-lattice relaxation with a 
related relaxation parameter Fi = T^^. And to take into account the spin-spin attenuation, 
we include the corresponding transverse width F2 = T2^- 

In this way, for the variables (43) to (45), we obtain the evolution equations 

— = -? (cjo + ^0 - i^2) u + fs , — = -2T2W + {u* f + f* u)s , 

f = -^(^*/ + r^)-ri(s-c), (47) 

where ujq is the Zeeman frequency (34) and C ^ [^Ij 1] is an average spin polarization. In 
the absence of pumping, ( = —1. But if there is a kind of pumping, for instance, by means 
of dynamic nuclear polarization, then ( > —1. 

Equations (47) are stochastic differential equations, containing the stochastic variables 
^0 and ^. The presence of the latter will make it possible to consider quantum effects. Also, 
Eqs. (47) are nonlinear due to the resonator feedback field entering through the effective 
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force (22). The resonator field is defined by the integral presentation (20), in which the 
magnetization density (13) writes 

"^x = - P /io Hu* + u) (^p = — j . (48) 

Thus, the spin evolution is characterized by the system of stochastic nonlinear integro- 
differential equations (20), (22), (47), and (48). These equations are assumed to be compli- 
mented by the initial conditions ■u(O) = Uq, w{0) = Wq, and s(0) = Sq. 

To proceed further, we notice that there exist several small parameters. Thus, the spin- 
lattice and spin-spin attenuations are usually small as compared to the Zeeman frequency, 
similarly to the dynamic broadening, 

— <1, — <1, — <1. (49) 

Uq Uq UJq 

The last inequality means that the infiuence of local random fields can be treated as weak, 
since the stochastic averages (41) are proportional to T^. The interaction energy of magnetic 
moments /i„ = h'jnl with the resonator field is proportional to the natural width 

^o = T VPfJ'P = T^VPlnl^n , (50) 

h 

which, being of order r2, is also small with respect to Uq. And the resonant circuit is assumed 
to be of good quality, i.e. Q ^ 1, because of which the resonator ringing width 7 is much 
smaller than the circuit natural frequency 00. Hence, there are two other small parameters 

■p 

— <1, -<1. (51) 

Specifying the external transverse magnetic field, entering the total field (3), we take 

Bi = ho + hi cos out , (52) 

where ho = const. And let the resonant part of the reduced electromotive force (17) be 

h = /i2 cos cut . (53) 

Defining the quantities 

Poho Pohi Afo^2 /-^N 

'^°"^' "^"^' '^^"^' (''^ 

we keep on mind that the amplitudes of all transverse fields (52) as well as (53) are small, 

in the sense that 

I t'n I I ^1 1 I ^2 1 / N 

^^<1, ^<1, ^<1. (55) 

UJO UJo ^0 

The existence of the small parameters (49), (51), and (55) allows us to realize the clas- 
sification of relative quasi-invariants. From Eqs. (47) we infer that the transition function 
(43) has to be considered as fastly varying in time, as compared to the slow functions (44) 
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and (45). This implies that w and s are temporal quasi-invariants with respect to u. Fast 
variation means that with a frequency of order ujq or u. These two frequencies are supposed 
to be close to each other, satisfying the resonance condition 



A 
-< 1 

LiJ 



U -UJQ , 



(56) 



which can be always done by tuning the natural resonator frequency to the Zeeman frequency. 
The occurrence of the small parameters, listed above, allows us to solve the integral 
resonator equation (20) by an iteration procedure, starting with the solution of the first of 
Eqs. (47) taken in zero order with respect to small parameters, that is, substituting in the 
right-hand integral of Eq. (20) the form u ~ ■Uoexp(— ic^^ot). Accomplishing the integration 
gives the first-order approximant 



h 



i{au — a*u*) + 2/?coscjt 



in which a = a{t) is the coupling function of spins with the resonator feedback field. 



a 



Toc^o 



1 



,-jAt-7t 



7 + iA 
and (3 = f3{t) is the coupling function of spins with the electromotive force. 

The resonance condition (56) is used in deriving Eqs. (57) with (58). 
An important quantity that appears is the spin-resonator coupling 



(57) 



(5^ 



(59) 



9 



7roc^o 



12(72 + A2) ■ 
This parameter enters the real and imaginary parts of the coupling function (58), so that 



(60) 



Re q; = gT2 
Im q; = —gT2 



sin At H cos At e^'^* 



A 

.7 V 7 

These formulas can be simplified if the resonance is sharp, such that 



cos At sin At J e~^* 



A 



|A| lAI 

7 U! 

Then the coupling function (58) becomes 

a = gT2 (1 - e--" 
This simplification is not principal but just helps us to avoid cumbersome expressions. 



(61) 



(62) 
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Taking into account in the effective force (22) the expression (52) for the transverse 
magnetic field and Eq. (57) for the resonator field, we define 

/i = —iuo — 2i{ui + l3) cosut + ^ . (63) 

Then Eqs. (47) can be transformed to the system of equations 

— = -i{uJo + ^o)u - {^2 - Ois)u + fis - asu* , 

— — = — 2(r2 — as)w + (m*/i + flu)s — 2Q;sRe v? , 

r/s 1 

- = -aw- - {u*h + flu) - ri(s - C) + aRe u^ . (64) 

Recall that, according to the classification of relative quasi-invariants, the function u is 
considered as fast, being compared with the functions w and s. The latter are temporal 
quasi-invariants with respect to u. The time derivatives of the coupling functions (59) and 
(62) are proportional to 7, because of which a and (3 are also quasi-invariants with respect 
to u. 

Following further the scale separation approach [9,42,56,57], we may solve the system 
of equations (64) by the method of stochastic averaging^ which is a generalization of the 
averaging technique [43] to stochastic and partial differential equations. First, we solve the 
equation for the fast variable u from Eqs. (64), with all quasi-invariants kept fixed. This 
yields 

I exp <^ -{iujQ + T2- as)t - i / ^o(^') dt' \ + 

+s I /i(t) exp {-{iujQ + V2- as){t - t') - i f ^o{t") dt"\ dt' . (65) 

Then, we substitute the solution (65) into the equations for the slow functions w and s, 
which are the second and third of Eqs. (64). After this, we average the resulting equations 
over time in the infinite interval, keeping the quasi-invariants fixed, and over the stochastic 
fields, employing the averages (41). To slightly simplify the resulting equations, we take the 
initial condition for the transition function u in the real form 

uo = j < /r(0) > . (66) 

Note that the counter-rotating term, containing u*, in the first of Eqs. (64) gives a small 
correction to the solution (65), of order r2/cJo- Hence, this negligible correction is omitted. 
In accomplishing the averaging of the second and third of Eqs. (64), the nonresonant terms, 
proportional to Re u"^, also give small contributions to the right-hand sides of these averaged 
equations. Taking all this into account, we could omit in Eqs. (64) the last terms and also 
could take the force (63) in the form 

/i = -zz/o-^(^^i+/9)e-*-* + e. 

This is a kind of the rotating-wave approximation or, in other words, the resonance approx- 
imation. 
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U = Uqi 



Recall that, when averaging the right-hand sides of Eqs. (64), we keep all quasi-invariants 
fixed. The function exp(— Ft) is also considered as a quasi-invariant. The stochastic average 
of Eq. (65) is 



<< u » = - 



UqS 



uJo — iT 
where the collective width 



[ui + p)s ^_^^^ 



A + ir 



Mo 



UqS 



ujq — iV 



A + ir 



-(iaJo+r)t 



T = V2 + T^- a s 
is introduced. Let us also define the effective attenuation 

P^P I ^gr _ Mvi + m -Tt I {^i + /^)^r . 



-,-rt 



(67) 



(6^ 



(69) 



Following the described procedure of averaging the second and third of Eqs. (64), we 
come to the equations for the guiding centers, 



dw ,^ , ~ o ds ~ ^ , ^, 

— = -2(F2 - as)w + 2F s^ , — = -aw -Ts- Fi(s - C) 

CLv (jjL 



(70) 



These are the main evolution equations describing the nonlinear dynamics of nuclear spin 
motion. The exponential factors entering the effective attenuation (69) characterize retarda- 
tion effects that are always present in real processes and, hence, are to be taken into account. 
This retardation may essentially influence the spin dynamics. 

6 Incoherent Quantum Stage 

One of the most interesting questions in the theory of spin superradiance is: What initiates 
the motion of spins when no transition coherence is imposed on the system at t = and no 
external flelds push spins in the transverse direction? That is, what is the origin of pure spin 
superradiance? In other words, how the coherence in spin motion develops from initially 
uncorrelated spin fluctuations? 

Let us note that here we are interested in the transition coherence that is related to the 
arising coherent radiation. In general, one may distinguish two types of coherence, state 
coherence and transition coherence [60]. A spin ensemble possesses state coherence, when it 
is prepared in a spin polarized state, with a nonzero spin polarization (45). And transition 
coherence develops when the transition function (43) is nonzero, or the coherence intensity 
(44) becomes noticeable. These functions (43) to (45) may be considered as dynamical order 
parameters [97] characterizing the level of state coherence or transition coherence. If the 
spin sample is initially polarized, it possesses state coherence. But if no external flelds are 
thrust upon the system, there is no transition coherence. An intriguing question is: How 
the transition coherence develops from the state coherence in a self-organized way? 

It is well known how transition coherence appears in a system of inverted atoms [4,5]. 
The relaxation begins with atomic spontaneous radiation, which is a quantum incoherent 
process. After the appearance of the seed radiation fleld, atomic correlations start arising 
through the interatomic photon exchange. Transition coherence develops as soon as the 
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interatomic correlations become sufficiently intensive. Then the quantum stage of spon- 
taneous emission changes for the coherent stage, resulting in atomic superradiance. Since 
moving spins produce magnetodipole radiation, one may ask if the latter can be the cause 
of collective spin motion, by analogy with what happens in atomic systems. 

6.1 Radiation of Magnetic Dipoles 

The vector potential created by radiating objects is 



Arad{r,t) = - J 3 fr',t 



r — r 



dr' 



c / I r — r' I 

where the density of current, formed by the magnetization density, writes 

j(r,t) = cV X m(r, t) . 



(71) 



(72) 



Let this current be produced by N^ff spins in the volume -^e//- These spins can act as an 
altogether only if the radiation wavelength is much larger than Le/f, so that 



koLeff < 1 



kn 



(73) 



The self-action of a radiating spin corresponds to N^ff = 1, with Lg// being the nucleus 
radius. 

Under condition (73), the vector potential (71) can be presented as an expansion 



Arad{r,t) ^- (- 

c J \x 



Id X d"^ 



X 



Q3 



c dt^2c^dt^ Qc^ 5i3|J(r + ^'^)^^ 



(74) 



in which x = r' — r and x = |x|. Assuming that there is no current through the surface of 
the volume L^j-j, for the magnetic field Urad = V x Arad, we have 



iirad{r, t) 



x-^ 



1-92 1 93\ ., 

+ TT^ 777^ I |X X j(r + X, t)J dX . 



2xc^ df^ 2c3 dt^^ 
The latter, involving Eq. (72) and defining the total magnetization 

M = — r X \ dr = /mrfr, 
2c J J 

can be transformed to 



(75) 



Hrad(r, t) 



3(m ■ x)x 



x^ 



m 



x"" 



(ix 



2c2 



rnxx m 

■^ \ 

x^ X 



(ix 



2 d^M 



3c3 dt^ 



(76) 



(77) 



where the overdot implies the time differentiation and m = m(r + x, t). 

Note that this way of presenting the magnetic field created by a system of radiating 
magnetic dipoles goes back to Ginzburg [98] and since then has been considered by many 
authors (see e.g. [99,100]). The first term in Eq. (77) describes a dipolar demagnetizing 
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field, depending on the shape of the sample. For the spherical shape, or for a sufficiently 
large volume, this term is zero, similarly to the first of Eqs. (6). The second integral in Eq. 
(77) can be approximately presented [100] as 



[m ■ xjx m 

x^ X 



8 (PM. 

dx. = 

SirLeff df^ 



Then one has 

4 •■ 2 d^M. 

Snc'^Leff 3c^ dt^ 

Here the term with three time derivatives corresponds to the so-called electromagnetic fric- 
tion. Generally, such terms with odd number of time derivatives lead to some kind of 
relaxation in the spin motion. For example, the term M in the Landau-Lifshitz equation 
[101] can be connected [102] with spin relaxation due to spin-phonon interactions. 
The radiation field (78), in which 



N, 



eff 



M = /io ^ < I. > , (79) 



1=1 



has to be added to the total magnetic field (3). Then in the Zeeman term (2) there appears 
the additional interaction — /ioH^ad-Ij- Wishing to concentrate the consideration on the role 
of the radiation field (78), we assume that there are no transverse fields, that is, -Bi = and 
H = 0. Then the Zeeman interaction (7) can be written as 

Hi = -/io {Bo + H^^^) I- - - fiQ (iJ+^ Ir + /+ H~^ , (80) 



where H^^^ = H^^^±iH^^^ is expressed through Eq. (78) and the magnetization components 

in which M~^ is the complex conjugate of M~ and /i„ = ^qI = h'-fnl- Introducing the 
notation for the radiation width 



2 
3 

we may write 



^rad = - k^'^nf'nNeff , 



1 „_ _ Trad (d\ 2C ..\ 1 _ Trad ( dh 2c 



Jeff ) It. UJQ \UL" '^J^eff 



Now, the evolution equations for the functions (43) to (45) can again be presented in the 
form (47), in which ujq is to be replaced by uq + fioH^^^/h and 



f — ~ T f'oHrad 



I 

n 

Keeping in mind the existence of the small parameter 

J- rad 



UJQ 
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« 1 , (82) 



we can find the radiation field (78) by iterating its right-hand side with the zero-order 
approximation u ~ Uq exp{—iujQt) and s ~ Sq. Then 

/ = Trad 1 r~F " 

and H^ad = 0. As a resuh, instead of Eqs. (47), we obtain 

du , 

— = -Z[UJQ + dujQ + 4o - lT2>U + TradSU , 

dw ,^ ^ , ds 

— = -2{T2-TradS)W , — = -TradW - Ti{s - Q , (83) 

with the frequency shift 

5UJQ = 



■^-L radS 



nkoLeff 



Equations (83) show that the magnetodipole radiation can lead to the arising collective 
effects only if the radiation width (81) is much larger than the transverse dephasing width 
r2 due to dipole-dipole interactions. For the latter one has 

r2 = riop-fnl^n , (84) 

where no is of the order of the nearest-neighbour number [12,23,53]. Comparing Eqs. (81) 
and (84), with taking account oi p = Neff/Ll^^, we find 

^™' ^ (A;oL,^/)^ « 1 , (85) 



Ta 3no 

since no ~ 10 and k^L^ff ^ 1 according to condition (73). Hence Trad ^ r2 and, respec- 
tively, the radiation time Trad = ^rad i^ much larger then the dephasing time T2 = F^^. 

To estimate the related parameters, we may take the typical value uoq ~ 10^ Hz, so that 
k^ ~ 10~^ cm~^. For protons, the gyromagnetic ratio 7„ = 2.675 x 10^ G~^s~^ and spin 
I = 1/2. The proton magnetic moment /i„ = 2.793/iAr can be expressed through the nuclear 
magneton /xtv = 5.051 x 10^^^ erg G^^. We have /i„ = 1.411 x 10^^^ erg G^^. The dimension 
of Gauss is such that G^=erg-cm^^. In this way. Trad ~ lO^^^^e// s^^, while F2 ~ 10^ s^^ for 
p ~ 10^^ cm~^. The relaxation caused by the self-action through magnetodipole radiation, 
when Neff = 1, yields the radiation time Trad ~ 10^^ s. This, with 1 year ~ 10^ s, gives 
Trad ~ 10^^ years. Such an enormous time is not only much larger than T2 ~ 10^^ s, but also 
surpasses the Earth lifetime of 5 x 10^ years and is even longer than the Universe lifetime 
of 10^*^ years. Clearly, the self-action, caused by the magnetodipole radiation, is completely 
negligible. Even with the account of collective effects, when N^ff ~ 10^^, we get Trad ~ 10~^ 
s~^, hence Trad ~ 10^ s, which is much larger then T2, their ratio being Trad/T2 ~ 10^. These 
estimates show that the ratio (85) is extremely small, Trad/T2 < 10~^. 

Thus, we come to the conclusion that magnetodipole radiation is absolutely unable to 
organize the coherent motion of spins. This conclusion is based on the inequality (85), which 
is always valid, independently of the nature of the spins involved, whether these are nuclear, 
electron, atomic, or molecular spins. 
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6.2 Resonator Nyquist Noise 

For many years, practically all researches have used to write that the major cause producing 
spin motion, resulting in pure superradiance is the thermal Nyquist noise of the resonant 
electric circuit. This belief has been especially surprising because Bloembergen and Pound 
[7] mentioned that this noise cannot be a noticeable cause of spin relaxation. The common 
delusion about the principal role of the Nyquist noise was clarified in the papers [10,55-57]. 
The role of this noise can be studied by analysing the effective attenuation (69). Assume 
that there are no transverse external fields, so that z/q = ^i = 0, but there is only the electro- 
motive force corresponding to the resonator Nyquist noise. Then the effective attenuation 
(69) is 

f = Ts + Tres , (86) 

where, in accordance with Eq. (59), 

The latter attenuation is caused by the Nyquist noise of the resonant electric circuit. For 
concreteness, we shall keep in mind the standard situation, when the resonator width 7 
is larger or of the order of r2 and T^. And, for simplicity, we consider the case of good 
resonance, when condition (61) is valid. Then at short time, the attenuation (87) reads 

Tres ^ \ yh2t^ {it < 1) . (88) 

Notice that Vj.es = at t = 0. 

The resonator electromotive force, entering the Kirchhoff equation (11), can be written 
as 

Ef = Eres cos out . (89) 

The amplitude of the reduced electromotive force (17) or (53) is /12, for which we have 

cEres ,2 ^^VEres 

/i2 = — -. — , ri2 ~ 



The amplitude squared of the electromotive force (89) caused by the Nyquist noise, can be 
presented [18] as 

with ut being the thermal frequency. At radiofrequencies, one has u <^ ut- Then Eq. (90) 
can be slightly simplified to 

i?L--7i?^T (-«l)- (91) 

71 \UJt / 



For the quantity z/2, defined in Eq. (54), we get 

2 _ 2TQEres ^ 2ToiJT 
^^ ~ h-fIRN ~ TilN 
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(92) 



For the attenuation (88) at t ^ I/7, we find 

'■^^ - 47 " 27r/7iV • ^''''^ 

This is, actually, the maximal value of the relaxation that can be reached by the attenuation 
(87). 

To estimate Eq. (93), we again accept the typical values of r2 ~ Fq ~ 10^ s ^ and 7 ~ 10^ 
s~^. At temperature T ~ 0.1 K, one has ujt ~ lO"*^^ Hz. Thus, we obtain F^es ~ (10"'^''/A^) s~^. 
If iV ~ 10^^, then T^es ~ 10~^^ s~^, which is much smaller than F3 ~ F2. The related time, 
during which the Nyquist noise could produce spin relaxation, Tres = T^es; is Tres ~ 10^^ 
s ~ 10^ years. This is many orders longer that T2 ~ T3 = F3"^. Therefore, the resonator 
Nyquist noise plays no role in spin relaxation of macroscopic samples and can never serve as 
a triggering cause starting the spin motion, unless the number of spins is much smaller than 

6.3 Local Spin Fluctuations 

We continue considering the case when there are no transverse external fields pushing spins, 
so that z/q = z/i = 0. As we have found out, neither the magnetodipole radiation, nor the 
resonator Nyquist noise are able to start the spin motion. So, we set z/2 = 0, hence P = 0. 
Then the effective attenuation (69) becomes F = F3. Recall that the width F3 is due to the 
dynamic broadening caused by local spin fluctuations, which are a kind of local spin waves. 
These spin fluctuations are, thus, the sole possible origin that could trigger the spin motion, 
when there are no external fields. 

To understand how the spin motion starts, consider short times, when 7t ^ 1 and the 
resonator coupling function (62) is yet small, a ~ 0. Then Eqs. (70), under F = F3, are 

dw o ds , ^ , , 

— = -2T2W + 2Tss^ , - = -(Fi + F3)s + FiC . (94) 

Their solution at short times reads 

/ F \ F 

wc^ iwo- Y' ^0 ) exp(-2F2t) + ^ Sq > 

This shows that, even if no transverse coherence is imposed on the spin system at the initial 
time, so that Wq = 0, the coherence, anyway, begins to arise triggered by the local spin 
fluctuations in the presence of a nonzero initial polarization. 

The resonator coupling function (62) increases with time. The incoherent quantum stage 
of spin motion lasts till that time when the coupling function a{t) reaches a value such that 
as = F2. Then, as is seen from Eqs. (70), fast generation of transverse coherence starts in 
the system. The crossover time tc, separating the incoherent quantum stage and the coherent 
regime of motion, is defined by the equality 

a{Qs{Q = F2 . (96) 
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Assuming that the crossover time tc is the smallest among other characteristic relaxation 
times, such as r, Ti, and T2, we can simplify the solution (95) to the form 

w{U) ~ w;o + 2r3tcSo , s(te) ^ So + 2rit,C , (97) 

taken at t = tc- This assumption about the time tc being the shortest is necessary for the 
existence of superradiance. From definition (96), we find the crossover time 

tc = r\J^^]^ (98) 

where r is the resonator ringing time. In order that the quantum stage could certainly 
change for the coherent regime, the crossover time (98) must be positive and finite, which 
requires that 

gs^ > 1 (te > 0) , (99) 

which is, actually, the condition of maser generation. In the case of sufficiently strong 
coupling, Eq. (98) reduces to 

tc-— {gso > 1) , (100) 

which shows that the crossover time can really be made shorter than other characteristic 
times. Since the spin-resonator coupling (60) is positive, condition (99) also tells us that 
the initial spin polarization sq is to be positive, which implies that the spin system is to be 
prepared in a nonequilibrium state, being in an external field satisfying condition (46). 

7 Regimes of Coherent Radiation 

After the crossover time (98), the coupling function (62) quickly increases and reaches the 
value a ~ gr2- At this stage, the motion of spins becomes coherent. The main regimes of 
coherent spin radiation are described below. 

7.1 Transient Spin Superradiance 

The transient stage corresponds to times larger than the crossover time tc but essentially 
shorter than the longitudinal relaxation time Ti. Then in Eqs. (70) one may neglect the lon- 
gitudinal relaxation parameter Fi. Consider the case, when there are no external transverse 
fields, so that the effective attenuation (69) is F = F3. Here we take into account that the 
resonator noise plays no role in spin motion. Being interested in the situation of sufficiently 
strong spin-resonator coupling (7 ^ 1, we see that gT2 ^ F3, since F2 ~ F3. Thence, we 
may omit in Eqs. (70) the terms containing F3. Thus, at this transient stage, we consider 

the equations 

dw ^ , . ds ^ , , 

— = -2F2(1 - gs)w , - = -gT2W . (101) 

These equations possess an exact solution [10,55-57] that reads 

w = (^] ' sech^ (^-^ , s = -^ tanh (^-^ + - . (102) 

\9^2j \ ^P J 9^2 \ Tp J g 
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The integration parameters in this solution are obtained by considering as the initial condi- 
tion the values (97) at t = tc, which gives the delay time 



^In 



rp + r, 



r — r 



(103) 



where 

r; = r^ + {gT2)^ {wo + 2T3tcsl) , r, = r^igso - 1) , r^r^ = i . (i04) 

In order that the delay time be larger than the crossover time, but finite, that is, tc < to < oo, 
it should be that Fp > F^ > 0, which yields 

Wo + 2Tstcsl > , 9So>l. (105) 

At the delay time (103), the coherence intensity is maximal, 

w{to) = (so - - ) (1 + 2F3t,) , s{to) = - , (106) 

\ 9) 9 

after which the transition coherence fastly decays, 

2 
w ~ 4w(to) exp(-2F„t) , s ~ -Sq + - (t > to) . (107) 

9 

The effective time of the coherent pulse, keeping in mind that Fstc ^ 1, writes 

To 



^ ^stcigsof 



:io8) 



^{gso - 1)2 + g^wo I {gso - l)^ + g^wo_ 

The necessary condition for superradiance is Tp < T2, from which one has the inequality 

(gso-lf + g^wo>l. (109) 

However, condition (109) is not sufficient for superradiance and defines three possible regimes: 

Wo = , gso > 2 {pure superradiance) , 

g'^Wo > 1 — {gso — 1)^ , gso > 1 (triggered superradiance) , 

g'^wo > 1 , gso < 1 (collective induction) . (HO) 

Recall that superradiance is a self-organized process, which requires that the crossover time 
tc be positive, hence, according to Eq. (99), it should be that gso > 1. Collective induction 
is not a self-organized process, but it is induced by an initially imposed coherence intensity 

Wo- 

Thus, there exist two types of transient spin superradiance, pure superradiance that 
is developing as a completely self-organized process and triggered superradiance when the 
process is mainly self-organized but at the beginning it is slightly pushed by an externally 
imposed coherence. As follows from Eqs. (107), it is only in the regime when gso > 2 that 
the spin polarization becomes inverted after the superradiant pulse. This may happen in 
either the regime on pure or triggered superradiance. When g ^ 2, the polarization s after 
to is almost completely inverted to the value — so- This effect of polarization reversal can be 
used for fast repolarization of solid-state targets employed in scattering experiments [29,103]. 
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7.2 Pulsing Spin Superradiance 

If there are no external fields acting on the spin system, then, after the transient superradiant 
burst, occurring at the delay time to, the transition coherence dies out. In order to produce 
a series of superradiant pulses, it is necessary to involve an external action. There are two 
ways of organizing a regime with a series of superradiant bursts, which can be called pulsing 
superradiance and punctuated superradiance. 

To realize the regime of pulsing spin superradiance, the inversion of spin polarization has 
to be supported by a permanent pumping, which can be accomplished by means of dynamic 
nuclear polarization. In that case, the longitudinal relaxation parameter Fi plays the role 
of the pump rate F^ that can be much larger than the spin-lattice attenuation of 10~^ s~^. 
The pump rate FJ can be 10 s~^ or larger. In what follows, we shall continue writing Fi 
instead of FJ, keeping in mind that the value of Fi corresponds to the pump rate. This is 
done in order to avoid cumbersome notations. In the presence of pumping, the stationary 
polarization ( is the pump polarization, which can reach the value oi ( = 1. 

We consider the case when there are no external transverse fields, so that the effective 
attenuation (69) reduces to F = F3. At larger time, when '-ft ^ 1, the coupling function (62) 
acquires its maximal value a = gT2- Then Eqs. (70) write 

^ = -2F2(1 - gs)w + 2F3s2 , J = -gV^w - T^s - Fi(s - C) • (HI) 

To find the conditions when the pulsing regime could arise, we need to consider the long- 
time behaviour of the solutions to Eqs. (111). For this purpose, we define the stationary 
solutions to these equations and accomplish the Lyapunov stability analysis [60]. The ap- 
pearance of complex characteristic exponents, related to particular fixed points, would mean 
the existence of oscillations around the corresponding stationary solutions. 

If the pumping is not very strong, so that g( '^ —1, then the stable fixed points are 

For the associated characteristic exponents, we find 

J^ - -ri (i - J^] , Jr - -2r.(i + |,ci) - ^ f 1 - ^' 



\ 2\gC\T2j ' ^' '""" Fi V 2F2. 

These exponents, being real and negative, show that the fixed point (112) is a stable node. 
When the spin-resonator coupling is weak or the pumping is not strong, so that \g(\ <^ 1, 
then the stable stationary solution is 

2 



(113) 




Ti + F3 V Fi 
The characteristic exponents, if Fi + F3 7^ 2F2, are 



^ 2FiF2(Fi-2F2-F3) 

' ^^^ '+(ri + F3)(Fi-2F2 + F3)^^' 
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4rir2r3 

which classifies the stationary solution (113) as a stable node. But in the special case, when 

ri + r3 = 2r2, (114) 

the characteristic exponents become 

Jt ~ -212 ± I i^V^V^gCf + ^ T,gC . 

Then solutions (113) describe a stable focus, if gC, > 0. This means that there exists a series 
of pulses approximately separated from each other by the separation time 



/2TiT; 



T.ep = vrW— ^ (0<(?C«1), (115) 

where T3 = T^^ and condition (114) is valid. 

For strong spin-resonator coupling and sufficient pumping, such that g( ^ 1, the sta- 
tionary solution to Eqs. (Ill) is 

116) 

The characteristic exponents 



W2 ^ 


Cri 


4c^ 


1 
9 ' 


4-- 


^^' 


+ r3)- 


T2T3 

qCFi 



show that solution (116) is a stable node. 

Now let us analyse the behaviour of stationary solutions to Eqs. (Ill) for arbitrary values 
of g(, but for different relations between the pumping rate Fi and the dynamic broadening 
F3. Thus, when Fi ^ F3, then 



* C Ti * > Ti 

J- 2J- 3 J^ 3 



i-(i + ^C)^ 

i 3 



;ii7) 



The corresponding characteristic exponents, if F3 7^ 2F2, are 

,, ^ 2(2F2 + F3)FiF2 , , ^ (4oC + 2F2 - F3)Fi 

J+ ~ -2F2 \ ■J^ 1 ^ gA J- ^ _Y _ K_ys_^ 2 3^^ 

' (2F2-F3)F3 ^^' ' 2F2-F3 

which defines the fixed point (117) as a stable node. And for the special case, when 

F3 = 2F2 , (118) 

one gets 



J± ~ -F3 ± 2 ^2gCT,Ts + 2 (^C - l)ri . 
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In such a case, the fixed point (117) becomes a stable focus, provided that gC, > 0. This 
imphes the existence of a series of pulses, with the separation intervals 



-^ Sep ^ 



/T1T2 



(Fi < Tg = 212) . 



;ii9) 



Finally, for a large pump rate, such that Fi ^ F3, we have the stationary solutions, 
which acquire different forms for different values of g(^. In the case of weak coupling, when 
gC < 1, we get 

F. 

Si ~C 1 



C^rs 



with the characteristic exponents 



;i-^C)ri 



(120) 



j+^-Fi, jf^ -2(1 -(?c)r2 

telling that the fixed point (120) is a stable node. For the intermediate case, when g( = 1, 
we find 



w^ ~ 



Fo 



/ 2VrJ' 






This point is again a stable node since 



j+ ^ -Fi + 2F2 M , Ji ^ -4r2 M 

V J- 1 V J- 1 

And for strong coupling and effective pumping, when g( > 1, we obtain 



W2 ~ 



(^C-i)ri 



I 



(K-2)F3 



(K-i)ri 



The related characteristic exponents 

J^ Fi 

J2 - 2^1 



1±^1-8^((?C-1 



are real in the interval 1 < g( < 1 + F1/8F2, but become complex if 



K>1 + 



8F2 



Under this condition, the frequency 



UJr 



is real. Then we may write 



"\ 


2T,T,igC 


-1- 


ST,/ 


Jt 


.-ir. 


±iuJoo 





(121) 



(122) 



(123) 



(124) 
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This means that the fixed point (122) is a stable focus. Here again there arises a series of 
pulses separated by the time 

(125) 

where, to simplify the expression, we take into account that Fi < r2. 

In this way, the regime of pulsing superradiance can appear under the conditions resulting 
in the formation of a series of pulses, which at long times, are separated from each other 
by the corresponding separation time (115), or (119), or (125), depending on the accepted 
conditions. This conclusion has been confirmed by direct numerical solution of Eqs. (Ill), 
displaying the regime of pulsing superradiance [9,63,65]. Experimental observation of this 
regime was done in Refs. [30-33]. Pulsing spin superradiance can be used for organizing 
superradiant operation of spin masers [60,65]. 

There exists another possibility of producing a series of superradiant pulses, which was 
named punctuated spin superradiance [66]. This can be done in the following way. Suppose 
that the transient spin superradiance has been realized, resulting in a superradiant burst 
at the delay time to. Then, according to Eqs. (107), at some time after to + Tp the spin 
polarization is reversed. For large coupling g, this reversal is practically complete. Now 
assume that we again inverse the spin polarization from — sq to sq, thus obtaining an inverted 
nonequilibrium system. Such an inversion can be realized in three ways: reversing the 
external magnetic field Bq, acting on spins by a resonant 7r-pulse, or turning the sample 180° 
about the x-axis. After the time to, counted from the moment when the newly nonequilibrium 
state is prepared, another superradiant burst will arise. When the second burst dies out, 
one can again invert the spin polarization by one of the described methods. Then one more 
superradiant pulse will occur. This procedure can be repeated as many times as required for 
creating a desired number of sharp superradiant pulses. Contrary to the regime of pulsing 
superradiance, analysed above, when the intervals between superradiant bursts are defined 
by the system parameters, the time intervals between superradiant pulses in punctuated 
superradiance can be regulated. Is is admissible to form various groups of pulses, with 
different intervals between separated groups. Hence, it is feasible to compose a code, like 
the Morse alphabet, which could be employed for processing information [66]. 

7.3 Induced Coherent Emission 

In the investigation of radiation regimes above, we have assumed that there are no external 
transverse fields that would permanently act on spins. The possible existence of such external 
fields has been reduced only to the preparation of initial conditions for wq and sq, or to 
supporting the value of the pumping parameter (. But the external transverse field (52) 
has been switched off. This was done for studying the process of self-organization, which 
would not be perturbed by external fields pushing spins in the transverse direction and, 
thus, helping to develop the transition coherence. The temporal arising of self-organized 
coherent motion of spins, from an initially incoherent state, is the most interesting and the 
least studied problem. Actually, the gradual appearance of the transition spin coherence has 
never been described before Refs. [8,60,66]. This is why we, first, paid the main attention 
to this problem. 
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Switching on the external transverse field (52) of course infiuences spin dynamics. The 
presence of this field makes the classification of possible regimes of coherent emission less 
evident. If the transverse field (52) is very weak, so that the leading part in the effective 
attenuation (69) is due to the dynamic dipolar broadening T^, then the analysed above 
regimes do not change much. But if the transverse field is sufficiently strong, so that in 
the effective attenuation (69) the term T^ is not the largest, but other terms are either also 
essential or even prevail over Fa, than spin dynamics may be noticeably disturbed. However, 
in general, one may distinguish two principally different regimes related to either weak or 
strong spin- resonator coupling. If this coupling is weak, such that qsq < 1, then even very 
strong transverse fields can induce coherent radiation with only Tp > T2. But when both 
the spin-resonator coupling is strong, so that qsq > 1, and the external transverse field is 
strong, then the pulse of induced coherent emission can become short, with Tp < T2. In any 
case, when the spin system is subject to the action of the transverse coherent field (52), but 
there is no pumping supporting the value of the pumping parameter ( > —1, then there can 
arise the sole main coherent burst, which may be accompanied by a series of small fastly 
attenuating oscillations. That is, the transverse field (52) induces only transient coherence 
radiation. This type of coherent radiation, essentially caused and infiuenced by an external 
transverse field, can be called induced coherent emission. This should not be confused with 
superradiance that is a self-organized process occurring without the action of any transverse 
field of type (52). 

As an example, let us consider how a strong constant transverse field would infiuence 
spin radiation. When in the transverse field (52) only the constant part Hq is present, but 
hi = 0, then the effective attenuation (69) takes the form 



t'n 



LJn 



Y2 



with uq defined in Eq. (54) and F in Eqs. (68) and (62). From here, we see, first, that the 
applied constant field has to be rather strong in order to lead to a noticeable, as compared 
to F3, term. Generally, the second term here is smaller than F3, if condition (55) holds true. 
If so, then the following process is not much different from what has been analysed above. 
Just in all formulas, we need to shift F3 replacing it by 

p . Z/o(F2 + F3) ^ Uq 

Such a shift would mainly influence the value of the delay time (103). In particular, the 
natural magnetic anisotropy of spin samples is often modelled by an external constant fleld. 
Hence, the existence of such an anisotropy would also play its role in triggering the initial 
spin motion [61,62]. However, the term F3 plays the major role. 

In the case, when the transverse fleld (52) consists of only a resonant alternating fleld, 
with the amplitude hi, then the effective attenuation (69) becomes 



where ui is given in Eq. (54), F = F2 + F3 — as, and the coupling function a being deflned in 
Eq. (62). We have thoroughly analysed the behaviour of solutions to Eqs. (70) for various 
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system parameters, when the resonant transverse field is present. This has been accomphshed 
by numerically solving Eqs. (70) for different initial conditions Wq and Sq. Note that the term 
in r, due to the resonant field, is of the order of z^f/r2, which can be comparable with or even 
larger than T^. Thence the resonant field infiuences spin motion essentially stronger than 
a constant transverse field, which is, actually, rather understandable. Numerically solving 
Eqs. (70), we found that, varying the system parameters and the initial conditions, it is 
possible to realize a variety of transient coherent pulses. Increasing the amplitude of the 
resonator field makes the duration of the coherent pulse longer. In all the cases the solutions 
w{t) and s(t), independently of their initial conditions, dies out to almost zero on the scale 
not longer than T2. The attenuation of the solutions slows down when increasing Ui. The 
latter, under conditions (55), is limited by the relation z/i/r2 ^10^. If the external resonant 
field is so strong that conditions (55) are no longer valid, then, instead of Eqs. (70), one 
should go back to Eqs. (47). In the latter case of a very strong resonant field, the solutions 
to Eqs. (47) may display chaotic behaviour [33,104,105]. 

In conclusion to this section, discussing different regimes of coherent spin radiation, it is 
useful to recall how the latter could be measured. An ensemble of coherently moving nuclear 
spins generates magnetodipole radiation with the total intensity 

/(t) = A |M(t)f , (126) 

in which 

N 

M = ^7„^ <I.(t)> 

is the total magnetization of N nuclei. The standardly considered quantity is the radiation 
intensity averaged over fast oscillations. For the intensity (126), this gives 

m = ^ filu;',NMt) , (127) 

where /i„ = h'-fnl and w{t) is the function of coherence intensity (44). The proportionality of 
the radiation intensity (127) to the number of spins squared is characteristic of any coherent 
radiation, which is not necessarily superradiance. The level of radiation intensity (127) is 
rather weak. Thus, the maximal intensity, when w{t) ~ 1, for proton spins, with /i„ ~ 10~^^ 
erg G-\ ujo ~ 10^ Hz, and A^ ~ lO^^, is only 7 ~ 10"^ W. 

However, despite such weak radiation intensity, it can be easily measured. This is because 
one can measure the power of the current 

Pit)=Rfit), (128) 

which is generated in a coil by the radiating spins. Employing the relation (14) between the 
electric current and the induced magnetic field, we get 

fit) - ^, HHt) . 

The resonator field H{t) can be found from Eq. (57). Neglecting the thermal Nyquist noise 
implies (3 = 0. Averaging over fast oscillations yields 



Yn H\t) = 2a\t)w{t) . 
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Using this, for the averaged current power (128), we have 

P{t) = gV2inuN (l - e-^*)^ w{t) . (129) 

Comparing Eqs. (129) and (127) for 7t ^ 1, we find 

Pit) 3QX^ 



lit) Stt^K 



(130) 



where A = 2t[c/uj. The ratio (130) can be quite large. For instance, under uj ~ 10^ Hz, hence 
A ~ 10^ cm, Q ~ 100, and Vc ~ 10 cm^, one gets the order of 10^. This explains why even 
a very low radiation intensity can be easily measured, as is mentioned in Section 2. 

8 Electron-Nuclear Hyperfine Coupling 

Since real condensed matter contains both nuclei as well as electrons, it is important to 
understand how the latter could infiuence the nuclear spin dynamics. Electron and nuclear 
spins interact with each other through hyperfine forces. A combined system of nuclear and 
electronic spins is described [53] by the Hamiltonian 

H = Hnuc + Hele + H^yp , (131) 

consisting of the parts related to nuclei, Hnuc, electrons, Heie, and their hyperfine interactions, 

IP- 

The nuclear Hamiltonian 



^hyp 



Hnuc = 2_/ -^« + 9 2_/ -^«i (132) 

is the same as in Eq. (1), with the Zeeman term (2) and the dipolar term (4). The electronic 
spin Hamiltonian is 

Hele = - 2 51 -JiJ ^i ' ^J + %e XI ^ " ^» ' (^33) 

in which Jij is an exchange interaction potential, Sj is an electron spin, and 7e is the electron 
gyromagnetic ratio defining the electron magnetic moment yUg = h'-feS = 0.928 x 10"^'^ erg 
G^^, which approximately equals the Bohr magneton hb = 0.927 x 10~^° erg G~^. The total 
magnetic field acting on electron spins is the same field (3) as that acting on nuclear spins. 
The dipolar interactions between electrons are much weaker than their exchange interactions. 
The Hamiltonian of hyperfine interactions is 

H,yp = ^El^ ■ S. + ^E E< ir S^ ■ (134) 

The first term here is the Fermi contact hyperfine interaction between nuclei and s-electrons, 
with the energy 

^ = Y^V7e|^(0)r, (135) 
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where ip{r) is the electron wave function. Using the estimate |-?/'(0)p ~ S/Airr^, where 
tb = 5.292 X 10~^ cm is the Bohr radius, we have 

2/l%„7e 

A ~ ^ . 



r 



B 



Since /io = ^7n, for protons, with 7„ = 2.675 x 10'' G^-'^s^^ we get A/h ~ 10^ s^^. The 
second term in the hyperfine Hamiltonian (134) describes the dipolar interactions, with the 
dipolar tensor 

< = -n' ^ (5., - 3<,. ng) . (136) 






Similarly to Eqs. (6), one has 



It is again convenient to pass to the ladder spin operators Sf = Sj ±iS^. Then the electron 
spin Hamiltonian (133) can be written as 

i i^j 

where in the Zeeman term, one has 

Hf^ = n^^BoS^ + ^ ^7e (5i + H) {st + S-) , (139) 

and the interaction term contains 

ifg^ = -J., {StS- + StS^) . (140) 

The hyperfine Hamiltonian (134) can be presented in the form 

i i^j 

Here the first term is due to contact interactions, with 

^^^ =^A {l+Sr + IrSl + 21^3^) . (142) 

And the second term corresponds to dipolar hyperfine interactions, with 



-"ii — '^ij h ^i 7 -'i ^i 7 -'i ^i + f^ijh ^i + Pij-^i ^i + 



+a,, {itS^ + ItSj) + 4 (j/^r + jr^l) , (143) 

where we use the notation 

a,, ^ A- , a., = i (a- - lA^ , /3,^ ^ 1 (a-^ - A^ - 2z Al]) . (144) 
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From the first of Eqs. (137), we liave 

i(.+i) j{+i) 3(.+i) 

Tlie resonator feedback field is, as early, defined by the Kirchhoff equation (18), but with 

m. = ^Y.iln<n> -le<S:>) . (145) 

i 

With this magnetization density, the integral presentation (20) is again valid. 

Writing down the Heisenberg equations of motion, we employ the commutation relations 
for the nuclear spin operators If and also for electron spins, 

[Oj , Oj \ = 20ijO^ , [Oj , Oj \ = ±0jjOj , [j j , Oj J = . 

Generalizing Eqs. (21), we define the local fields acting on nuclear spins 

e = 1 E (2&^.^; - \ «..^7 + 2c.,^; + A,^; - \ ^^sj + a,,si] , (we) 

and local fields acting on electron spins, 

^0 = :^ E [^^■3^ + 4^7 + "^^-^i ) + ^ ( -^o^- - E -^^.^l 

^ = T E (A,^; - \ «.,/7 + ^^.-^l) + \ \J^S^ - E J^3Si\ , (147) 

where Jo = Z]j(^j) Jy- Introduce the frequency operators 

fi„ = -7„5o + - 5f + ^0 , ^e = leBo + 4 ^f + v^o , (148) 

and the force operators 

fn = -tln{Bi + H)+t^Sr+^, f^ = ^^^(B, + H)+^^Ir + ^,. (149) 

h n 

Then the Heisenberg equations of motion can be presented in the form 



(it "^ ^•"" cit 2 

^ = -^^eS■ +Sth, ^ = -1 {feSr + Stfe) • (150) 
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In the stationary regime, Eqs. (150) define coupled nuclear-electron spin waves, whose 
description can be done similarly to Section 4. To accomplish this, we may define, by analogy 
with (24), 

sr^<sr> +6sr . (isi) 

Here, as in Eq. (25), we assume an ideal lattice, for which < 5*^ >=< 5*° >. Considering 
the case without transverse magnetic fields, when i?i = if = 0, we set < S^ >= 0, while 
< Sf >^ 0. This is analogous to setting < Jj >= 0, while < If >^ 0. Hence, Jj = 61^ 
and Si = SSi , because of which the local fields (146) and (147) correspond to random local 
spin fiuctuations, with ^o = ^(.o, C = ^C^ Vo = '^V^o, and ip = Sip. Then the second and fourth 
of Eqs. (150) give (5/f = and 6S^ = 0. To the remaining equations, we substitute the 
Fourier transformed nuclear spin operators (28) and 

^f = E si exp (Ttk ■ r^) . (152) 

k 

We invoke the Fourier transforms (29) and also introduce 

ak= Yl '^ij 6xp {-ik ■ Tij) , Pk= Yl f^ij *^^P (^^'^ ■ ^ij) ' 

Jk= Y Jio exp (-ik ■ Fij) . (153) 

Generalizing Eqs. (31), we use the notation 

^^, ^ -7„5o + H < ^^ > +^ < ^^ > ' A, = ^ < If > , 

Sk = leBo + j-iJk- Jo) < S] > +^ < /| > . (154) 

Then Eqs. (150) can be reduced to the equations 

^^ = /i.47-A,/+- 1(a- ^) </;> S,- ^ <i:> Si, 

'^ = ^^S^- l{^- f) <S'> lu- I <St> n, (155) 

which describe coupled nuclear-electron spin waves. The spectrum of the latter follows from 
Eqs. (155) yielding the fourth-order algebraic equation defining two positive branches of 
coupled spin-wave collective excitations. The coupling between the branches comes from 
the terms of Eqs. (155) containing the hyperfine parameters A, a^, and j3k- Without the 
latter, there would exist two separate types of collective excitations corresponding to nuclear 
spin waves with the spectrum (33) and to electron spin waves with the spectrum 6^. In the 
presence of the hyperfine coupling, one may conditionally distinguish one of the branches of 
the coupled collective excitations as the branch related to nuclear spin waves, while another 
branch can be ascribed to electron spin waves. The hyperfine coupling, of course, essentially 
modifies the spectra of the coupled nuclear-electron spin waves. In particular, nuclear spin 
waves become noticeably infiuenced by the indirect exchange through electron spins, when 
the latter possess a long-range magnetic order. The related exchange force is called the 
Suhl-Nakamura force [54], which can stabilize the nuclear spin waves. When nuclear spins 
are initially in a nonequilibrium state, nuclear-electron spin waves play the role of a trigger 
starting the motion of nuclear spins from the nonequilibrium to their equilibrium state. 
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9 Enhanced Nuclear Radiation 

To describe the nonequilibrium dynamics of coupled nuclear and electron spins, we again 
employ the scale separation approach [9,42,56,57,106]. To this end, the variables (146) and 
(147), characterizing local field fiuctuations, are treated as random fields. The stochastic 
averages for the random variables (146) are defined as in Eq. (41). However, the dynamic 
broadening T^, in the presence of hyperfine interactions, is renormalized according to the 
relation 

'- 3 — '- nn~^ '- ne J (156) 

where r„„ ^ PnlnfJ^n is the broadening due to dipolar nuclear interactions, with p„ being the 
density of nuclear spins, and r„e ~ y/PnPelnfJ'e is the broadening caused by the hyperfine 
nuclear-electron interactions, with pe being the density of electron spins. 
For the random fields (147), the stochastic averages can be defined as 

< ipo{t) > = < ip{t) > = , < ipo{t)ipoit') > = 2735(t - 1') , 

< MtMt') > = < ^{tMt') > = , < ^*it)ip{t') > = 2735(t - 1') , (157) 

with the dynamic broadening 73 given by the relation 

7| = rL + rL. (158) 

Here, Fee ~ Pelepe is the broadening caused by electron spin fiuctuations, while Fe„ ~ 
\/PePnlePn is the broadening due to electron-nuclear hyperfine interactions. Recall that 
Hn = filnl and /ie = ^leS are the nuclear and electron magnetic moments, respectively. 

To estimate the related broadening widths, we may assume that Pn ^ Pe ^ 10^^ cm^'^. 
Then, with 7„ = 2.675 x 10^ G-h-\ 7e = 1.759 x 10^ Q-^fi'^ , Pn = 1-411 x lO^^s erg G-\ 
and pe = 0.928 x lO'^o erg G-\ we find F„„ ~ 10^ s~\ F„e ~ 10^ s'\ Fee ~ 10^° s-\ and 
Fen ~ 10^ s~^. As is seen, hyperfine interactions may strongly infiuence local nuclear spin 
fiuctuations but are not important for electron spin motion. 

By analogy with Eqs. (43) to (45), we introduce the functional variables for electron 
spins: the transition function 



1 ^' 



X 



SNe ,., 



E < ^7 > ' (159) 



where A'^e is the number of electrons, the coherence intensity 



and the electron spin polarization 

1 



y = o..wl l^ E < ^^^7 > ' (160) 



«7^i 



A^. 



SN,. 



E < 5| > . (161) 



The averaging in Eqs. (159) to (161) is accomplished over the spin degrees of freedom, not 
involving the random variables (146) and (147). The stochastic averages for the latter are 
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defined in Eqs. (41) and (157). The random variables ascribed to nuclear and electron 
fluctuations are assumed to be uncorrelated with each other. 

Employing expressions (159) to (161), we may present the averages of the frequency 
operators (148) as 

where the effective nuclear and electron frequencies are 

UJn = -InBo + ^ ASz , UJe = Te^Q + ^ AI S . (163) 

n n 

For an external field i?o ~ 1 T= 10^ G, since 7„ ~ 10^ G~^s^^, 7e ~ 10^ G^^s^^, and 
A/h ~ 10^ s~^, we see that the hyperfine term in the nuclear frequency cj„ can become 
larger than the nuclear Zeeman frequency cjq = |7n-Bo| ~ 10^ Hz, provided that there exists 
a nonzero average electron magnetization, with z ^ 0. Hence, it is necessary to accurately 
take into account the possible hyperfine shift of the nuclear frequency. At the same time, 
the electron Zeeman frequency 7e-Bo ~ 10^^ Hz is much larger than the hyperfine frequency 
shift, so that Ue ~ 7e-Bo- 
For the averages 

fn = <fn> , fe = < fe> (164) 

of the force operators (149), we have 

/„ = -27n(5i + H) + ^ASX + ^, /e = 27e(5i + H) + I- AIu + if . (165) 

Here, the hyperfine terms play the role of additional fields acting on nuclear and electron 
spins. 

It is worth mentioning that the hyperfine shifts in the effective frequencies (163) are 
caused by the first-order hyperfine effect. Generally, there exists the second-order effect due 
to the Suhl-Nakamura forces. The second-order frequency shift appears as follows. Let us 
consider the evolution equations (150) for the electron spin operators, with the explicit form 
of the variables (147), without treating the latter as random. Then, the third of Eqs. (150) 
yields 



where 






UJp. 



if 1 _ 

Substituting this to the first of equations (150), we come to the conclusion that the effective 
nuclear frequency should have the form 

OJn = OJo + J <St> +As7V , (166) 

in which the Suhl-Nakamura frequency shift is 

^SN=Y.W^ <St> s. (167) 
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This is often called the dynamic frequency shift because of its dependence on s = s{t). 

However, the second-order frequency shift (167) is much smaller than the first-order 
hyperfine shift in Eq. (166), 



A 

Thus, Eq. (167) shows that Asm ~ noSpnPe{ft-lnleY /^e, provided that there is a long-range 
magnetic order of electron spins, i.e., < S^ >^ 0. Here hq ~ 10 is the number of nearest 
neighbors. For the typical values of parameters, this gives Asn ~ 10^ Hz, which is much 
smaller than A/h ~ 10^ Hz, hence hAsN/A ~ 10~^ is really a very small value. This makes 
it admissible to take account of only first-order hyperfine frequency shifts. 

Averaging Eqs. (150) over the spin degrees of freedom, we shall use the decoupling of 
nuclear and spin operators 

< /"S-f > = < If >< S-f > , (168) 

valid for all i and j. And, similarly to Eq. (42), we employ the decoupling 

< S^S^ > = <S^>< S^ > (i ^ j) (169) 

for different indices i and j. Recall that this type of decoupling has been called the stochastic 
mean-field approximation, since the averages < . . . > do not include the averaging over the 
stochastic degrees of freedom, thus, allowing for the consideration of quantum effects [8,9]. 
Finally, from Eqs. (150), we come to the stochastic differential equations describing the 
evolution of nuclear spins, 

— = -i{LJn + 6 - i^2)u + fnS , — = -2T2W + {u* fn + fnU) S , 

ds 1 

^ = -2(«*/n + /»-ri(s-c), (170) 

as well as the motion of electron spins, 

— = -i{uJe + (P0- il2)x + feZ , "- = -2721/ + {x* fe + f*x)z , 

f^=-\{x*f^ + f:x)-^,{z-a), (171) 

where 71 and 72 are the spin-lattice and spin-spin relaxation parameters for electronic spins, 
while a is an equilibrium value of < S^ > / S. Equations (170) and (171) are to be compli- 
mented by the Kirchhoff equation (18), with rrix given by Eq. (145), which reads 

1 1 

?nx = - PnPm {U* +U) - - PePe {x* + x) . (172) 

The evolution equations for nuclear and electron spins are very similar to each other. This 
means that we could consider electron spin superradiance on the absolutely same footing 
as the nuclear spin superradiance. The major difference is that hyperfine interactions can 
essentially infiuence the motion of nuclear spins, while these interactions do not disturb much 
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the electron spin motion. Therefore, we concentrate now our attention on the evolution of 
nuclear spins. 

The system of equations (170) plus (171) can again be solved involving the scale separa- 
tion approach [8,9,42,56,57]. For this purpose, we remember the existence of several small 
parameters defined in Eqs. (49), (51), (55), and also assume the validity of the resonance 
condition (56). In addition, we have as well the inequalities 

^^<1, P^<1, r^<i- (173) 



\U)e\ \Ule\ \^e. 

Also, we have to keep in mind that 

^ = ^ « 1 . (174) 

Being based on these inequalities, we may realize the following classification of relative 
quasi-invariants: The functions w and s are temporal quasi-invariants with respect to the 
functional variables m, x, y, and z\ the functions y and z are quasi-invariants with respect to 
x; and u is also a quasi-invariant with respect to x. 

To derive the expression for the resonator feedback field, we again use the integral form 
(20), iterating it with the approximate solutions for x and m, keeping in mind that, because 
of the inequalities (174), we have \oOn/oOf,\ <^ 1. We substitute in integral (20) the solutions 

Xq - - — u e"*'^^* + - — u 

huJe / huJe 

and u = uoexp{—iuJnt). This shows that electron spins oscillate with the frequencies Ue, u;„, 
and uJe + uJn- Calculating the transverse magnetization (172), we insert it into integral (20), 
replacing z by its average a. In that way, we come to the expression '~fnH = fioH/h having 
the same form (57), but with the coupling function 

a = gV^ (l - e-^*; 
instead of a, where the renormalized spin- resonator coupling is 



9 = 9{'^ 






Here g is defined in Eq. (60) and uOe — 7e-Bo. Taking into account condition (46), according 
to which 7„-Bo < 0, the renormalized spin-resonator coupling can be presented as 

«.<,(l + ^). (175) 

where ujq = |7„i?o|- Analogously to the previous procedures, we may study different regimes 
of coherent nuclear radiation, just replacing everywhere the spin resonator coupling (60) 
by its renormalized value (175) and by using expression (156) for the dynamic broadening 



F = /F^ -I- r^ 

3 y -*- nn ' ne' 
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The noticeable renormalization of the couphng (175) occurs when there exists a long- 
range magnetic order in electron spins, so that a 7^ 0. Then the coupling (175) can become an 
order larger than (7, that is, g/g ~ 10. This, respectively, leads to the enhancement of nuclear 
spin radiation, with the current power (129) becoming an order stronger. Simultaneously, 
the crossover time (98), or (100), and the pulse time (108) become an order shorter. Thus, 
in the case of pure spin superradiance, when wq = 0, and for gsQ ^ 1, we have 



Then the delay time (103) changes to 



r 

gso 


Tp^ 


T2 

^So 


les to 




gso 


In 

2^ So 


2 



which shows that to also becomes about an order shorter. 

It is important to stress that the enhancement effect can exist not only in ferromagnets 
or ferrimagnets, but also in paramagnets, provided that the external field is sufficiently large 
and the temperature is low. Thus, for a paramagnetic system of electrons, the average 
magnetization per spin, due to the external magnetic field Bq, is 

a = tanh 



For the typical values of parameters, with Bq ~ 10"^ G, we have /ig-Bo ~ 10"^^ erg, hence 
fieBo/h ~ 10^"*^ s^^. The thermal frequency ujt = ksT/h, at T ~ 0.1 K, is ujt ~ 10^^ s^^, 
which is of the same order as fieBo/h ~ Ue- Then a ~ tanh(l) fti 0.76. Therefore, the 
enhancement effect, under such conditions, exists for a paramagnet as well. 

10 Superradiance by Magnetic Molecules 

There exists an interesting class of composite objects behaving as giant nuclei with total spins 
ranging from 1/2 to rather large values. These are magnetic molecules [60,107-109], which 
can form crystalline materials termed molecular magnets. In these materials, all magnetic 
clusters possess the same shape, size, and orientation, because of which the inhomogeneous 
broadening caused by the system nonuniformity is very low. As examples of such magnetic 
molecules, we may mention the following, where in brackets the total ground-state spin is 
shown: 

Ke [V^5+AS6042(H20)] ■ 8H2O (S = 1/2) 

[(PhSi02)6Cu6(02SiPh)6] (5 = 3), 

[Mni20i2(CH3COO)i6(H20)4] ■ 2CH3COOH ■ 4H2O {S = 10) , 

[Fe802(OH)i2T6]'+ (S = 10) , 

Mn604Br4(Et2dbm)6 (S = 12) , 

[Cr(CNMnL)6] (0104)9 {S = 27/2) . 
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The first of tlie above molecules is briefly denoted as V15. Molecular magnets formed 
by these molecules have no magnetic anisotropy. The short-hand notation for the third 
molecule is Mni2. The corresponding molecular magnets possess a rather strong single-site 
anisotropy characterized by the parameter D ^ 0.967 x 10~^^ erg, which gives D/ks ~ 0.7 
K and D/h ^ 0.917 x 10^^ s^^. At temperatures below the blocking temperature Tg ~ 3 K, 
the magnetization of a molecular crystal formed of Mni2 is preserved during the relaxation 
time Ti ~ 10^ s. The size of each molecule Mni2 is about 10 A and the distance between 
the nearest neighbours is around 14 A. Hence the average density is p ~ 0.364 x 10^^ cm~^. 
The spin-spin relaxation parameter is r2 = nQp^s\l^o\S, where 

fJ'O = -gsfJ'B = -^75 , (176) 

Qs is the Lande factor; hb, Bohr magneton, and 75 is the gyromagnetic ratio of a molecule 
with spin S. This parameter r2 ~ 10^*^ s~^ is due to dipolar spin interactions. 

The fourth molecule has the abbreviation Feg. In its chemical formula, the letter T stands 
for the organic ligand triazacyclononane. The corresponding molecular magnet possesses the 
magnetic anisotropy described by the parameter D ^ 0.414 x 10~^^ erg. Then, D/ks ~ 0.3 
K and D/h ^ 0.392 x 10^^ s^^. Below the blocking temperature Tb ~ 1 K, the relaxation 
of magnetization to zero in a molecular magnet occurs during the relaxation time Ti ~ 10^ 
s. The average density is p ~ 0.4 x 10^^ cm~^. The interaction between molecules is also 
through dipolar forces, resulting in the spin-spin attenuation r2 ~ 10^° s~^. 

In the last molecule from the list above, the letter L in its chemical formula means a 
neutral pentadentate ligand. Its properties are close to those of the molecules Mni2 and Feg 
described above. 

In principle, at very low temperatures, of the order of hT2/kB ~ 0.1 K, a purely dipo- 
larly interacting molecular magnet can become magnetically ordered. The paramagnet- 
ferromagnet transition at T^ ~ 0.16 K was observed [110] in the dipolar magnet formed of 
the molecules of the fifth type from the list above, this kind of a molecule being abbreviated 
as Mng. The molecular magnet composed of the latter molecules has a week anisotropy, with 
D ^ 1.795 X 10"^^ erg. Thus, D/ke ^ 0.013 K and D/h ^ 1.7 x 10^ s'^ The spin-spin 
dipolar relaxation parameter is r2 ~ 10^° s^^. The spin-lattice relaxation time, at Bq ~ 10^ 
G and T ~ 0.1 K, is Ti ~ 10^^ s, so that Fi ~ 10^ s^^. Generally, Ti can be estimated from 
the formula 

T,-roexp(^^-^j . 

For the molecule Mng, one has tq ~ 3 x 10"^ s. 

The Hamiltonian of a molecular magnet, consisting of N magnetic molecules, each having 
spin S, has the form 

^ = E^^ + ^E^^.. (177) 

in which Hi is related to individual spins and Hij, to pair spin interactions. The individual 
term 

Hi = -poB ■ S, - D{S^f (178) 

includes the part characterizing the single-site magnetic anisotropy with the anisotropy pa- 
rameter D. Positive D > implies an easy- axis anisotropy, while D < means an easy-plane 
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anisotropy. The pair term in the Hamiltonian (177) corresponds to dipolar spin interactions 

H., = Y.CtfStSl (179) 

a/3 



with the dipolar tensor 



All notations in this sections are close to those of Section 3. The difference is that now we 
are considering the molecular spins 5*^ formed by electrons. The ground-state molecular spin 
S can be quite large. The principal distinction, as compared to Section 3, is the necessity of 
taking into account the magnetic anisotropy that can be rather strong. As early, the dipolar 
tensor enjoys the properties 

Similar to Eq. (9), we use the notation 

with Eq. (10) being valid. Employing the ladder spin operators Sf = Sf ± iSf, we may 
present the term (178) as 

H, = -^i,B,St - D{Stf - i /io(i?i + H){St + Sr) , (180) 

where the total magnetic field is assumed to be the same as in Eq. (3). And the pair term 
(179) takes the form 

Hij = ttij ( D j Sj — -^ ^i ^j j + 
In the Kirchhoff equation (18), the magnetization density is 

m. = y E < ^r > • (182) 

i 

The fluctuating local flelds are denoted as 

■Co = 7" E y^ij'^^j + ^*ij'^j + ^ij'^j 

e = r E (2%^i - \ <^^^SJ + 2h,Sl) . (183) 

Writing down the Heinsenberg equations for the spin operators, we shall need, in addition 
to the standard commutation relations, the relation 
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Thus, the evolution equations for the spin operators, with the notation 

f = ^Js{B^ + H)+^, (184) 

become 

^ = -^(75^0 + Co)S- + fSt + I J {S-S^ + S^S-) , 

"^^^ - I (rSr + Stf) . (185) 



dt 2 

Analogously to Section 4, we can show that there exist molecular spin waves. A subtle 
point here is the linearization of the last term in the first of Eqs. (185), which has to be 
done so that to take into account the absence of this term for S = 1/2. This can be achieved 
by invoking the linearization 

S-S; + S!Sr^(2-^)<Sr>S- (186) 

possessing the correct asymptotic behaviour for S = 1/2 as well as for S ^ oo (see discussion 
in Refs. [60,94]). Introduce the effective frequency 

UD = 75^0 - j(2- ^^<S:> . (187) 



From Eqs. (185), we obtain 



dt 
where 



-i/ifcS'fc +iXkS^ 



/ifc = ^ < Sf > +a;z, , Afc = ^ < Sf > . (189) 

Then the spectrum of molecular spin waves is described by Eq. (33). The condition for the 
spectrum to be positive, in the case of a cubic symmetry reads 



2pnlSuJo 



< 1 (uJo = -\fioBo\] . 



The external longitudinal magnetic field is directed so that, similarly to the inequality (46), 
we have 

lioBo < , 75^0 > . (190) 

Following Section 5, we define the transition function 

1 ^ 

« = ;5]7E<^r>, (191) 



the coherence intensity 



1 ^ 
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and the spin polarization 

E<^r>, (193) 



1 ^ 



where the angle brackets denote the averaging over the spin degrees of freedom, interpreting 
the local fields (183) as random variables with the stochastic averages (41). For the products 
of spin operators with different indices, we employ the stochastic mean-field approximation 
(42), while for the operators with coinciding indices, we use the decoupling 

< S-S^ + S^S- > = (^2- i) < S- X St > 

resulting from Eq. (186). 

The uniform molecular spin- wave frequency (187) becomes 

u;d=u;o-{2S-1)j-s, (194) 

with uq = 75-Bo- Then for the functions (191) to (193) we derive the same equations (47), 
except that the first of these equations appears as 

— = -liuD + ^0 - iT2)u + fs , (195) 

with uq replaced by ud- Note that uq is always positive, while ud, given by Eq. (194), can 
be both positive as well as negative. 

To deduce the expression for the resonator feedback field, we resort to the integral Eq. 
(20). In analogy with Eq. (57), we find 

— — = i{a£)U — a})U*) + 2j3cosujt , 

lb 

where, instead of Eq. (58), we have the coupling function 



ao = Eqc^d 



1 — exp{— ^(c^ — cj£))t — 7t} 1 — exp{i(cj + ujD)t — 7t} 



in which the natural width is 



(196) 



To = TT'qp'jslJ'S {fJ'S = ^IsS) . (197) 

An efficient coupling with the resonator occurs only when u ~ uj£,, if uj£, > 0, or when 
uj ~ —ujd, if ujd < 0. This implies the resonance condition 

\Ad\ , , 

<^ 1 , Ad = 00 — \ujd\ ■ 



CO 

If the resonance is sharp, such that \Ad\ < 7, then the coupling function (196) reduces to 

«D = 90^2 (l - e-^*) , (198) 
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where the spin-resonator couphng is 

In this way, introducing the effective force 

we come to the evolution equations 

— = -i{ujD + Co)m - (r2 - aDs)u + fis , "17 = -'^{^2 - aDs)w + (m*/i + A* «) s , 

ds 1 

= ^anw - - (w7i + flu) - T,{s - () . (200) 

at 2 

Here we have omitted the nonresonant terms, corresponding to the last terms of Eqs. (64), 
which do not contribute to the equations for the guiding centers in the averaging technique. 
Using this technique, we obtain the guiding-center equations of the same form as Eqs. (70), 
but with a, uq, and A replaced by an, ujd, and A^, respectively. The following description of 
possible regimes of coherent radiation by magnetic molecules is the same as that for nuclear 
spins. 

First of all, it is necessary to stress that the electromagnetic interaction of radiating 
magnetic dipoles through the common radiation field can never produce coherent radiation. 
This was explained in detail in Section 6.1, whose results are valid for arbitrary spins, whether 
these are nuclear or molecular spins. The impossibility of collectivizing the spin motion by 
means of the magnetodipole radiation is caused by the fact that the corresponding radiation 
time Trad ^ ^2 is many orders larger than the spin-spin dephasing time due to dipolar spin 
interactions. This is expressed by inequality (85). Therefore, the assumption [111] that 
molecular nanomagnets could exhibit superradiance caused by magnetodipole radiation is 
wrong. In addition, the regime considered in Ref. [Ill] requires a very strong transverse 
magnetic field of about 6T~ 10^ G. Such a strong transverse field suppresses the single- 
site anisotropy, whose effective field is Bo = (2S — l)D/h'~fs- The latter, e.g., for Mni2, 
with 7s ~ 10^ G^^s^^ and D/h ~ 10^^ s^^, is of the order of Bd ~ 10^ G. Any magnetic 
relaxation in the presence of such a strong transverse magnetic field has nothing to do with 
superradiance, which is a self-organized spontaneous radiation. When a transverse field 
induces spin relaxation, this is called spin induction. Even if a kind of spin induction could 
be realized for magnetic molecules, the related relaxation time would be of the order of T2. 
So, it would be just a simple spin induction, without any collective effects. 

In order to realize the real spin superradiance by magnetic molecules, it is necessary to 
couple them to a resonant electric circuit, whose feedback field could collectivize the spin 
motion. When there are no transverse external fields, as it should be for the regime of pure 
superradiance, the process is triggered by molecular spin waves. Then the arising feedback 
field can make the spin motion coherent, in the same way as it has been considered above 
for nuclear spins. The difference with the latter is the presence of magnetic anisotropy that 
may complicate the experimental realization of molecular spin superradiance. 
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From the evolution equations (200) it is seen that the increase of the transverse coherence 
occurs when the effective attenuation r2 — Qds becomes negative, which requires that the 
spin-resonator couphng (199) be positive and sufficiently large. To have positive qd, one 
needs that the effective frequency (194) be also positive, that is 

ujo> (2S-1) — s . 
n 

This happens, assuming that s > 0, when D is negative, so that the molecules form an 
easy-plane nanomagnet, or, if D > 0, when the applied longitudinal field Bq is sufficiently 
strong. For instance, in the case of Feg, the anisotropy field is Bd ~ 10^ G. Hence, the 
external field has to be stronger than 1 T. 

Another complication is that the frequency (194) is, actually, a function of time through 
s = s{t). Therefore, to organize the resonance condition uj ~ ujo, one has, in general, to 
vary either the resonator natural frequency uj or the external field Bq, so that to achieve the 
approximate equality 

^sBo = uj + i2S-l)^s iAD = 0). 

n 

In principle, such a temporal varying of the system parameters, for achieving the resonance 
conditions, is feasible and is known for optical systems [5], where it is called chirping. The 
chirping technique can also be used for realizing superradiance from magnetic molecules [60]. 

The situation becomes simple when S = 1/2, so that the anisotropy disappears, or when 
ujq ^ (25* — l)D/h. Then uj£, ^ ujq, and all the consideration reduces to the same as has 
been done above for nuclear spins. Another possibility is by tuning the resonant electric 
circuit to one of the transition frequencies of admissible 25* transitions and by supporting 
the population of the upper level with the help of the permanent pumping, as it was done 
for the nuclear spin / = 5/2 of ^^Al in experiments [30-33]. In the latter case, solely the 
regime of pulsing spin superradiance can be achieved. 

In this way, though the presence of magnetic anisotropy complicates the experimental 
realization of coherent spin radiation by magnetic molecules, nevertheless there are several 
possibilities for reaching the required conditions and producing such a coherent radiation. 
This could be used for superradiant operation of spin masers [60] . Radiation from magnetic 
molecules would be essentially stronger than that from nuclear spins. And not only the 
current power (129) could be easily measured, but radiation intensity (127) would also be 
very high. To estimate the radiation intensity (127), let us take ujd ~ cuq ~ 10^^ Hz, which 
corresponds to microwaves with the wavelength A ~ 0.1 cm. For the parameters typical for 
Mni2 or Feg, we have r2 = noP7s/is ~ 10^° s~^ and Fq = O.IF2. The spin-resonator coupling 
(199) is quite large, go ~ 10^. The radiation intensity / ~ ^^gu'^N'^/c', for A^ ~ 10^°, 
reaches the value I ~ 10^^ W. The superradiant pulse can be very short. The pulse time 
Tp ~ T2/ gso can be as short as 10^^^ s. So, it is feasible to realize femtosecond superradiant 
pulses. The possibility of producing coherent radiation by magnetic molecules can find 
various applications. 
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11 Pion and Dibaryon Radiation 

Nuclear systems under extreme conditions can emit not only photons but also other par- 
ticles. If the emitted particles are bosons, they can form coherent beams, similarly to the 
coherent photon radiation. The required extreme conditions can be experimentally achieved 
in the process of nuclear and heavy-ion collisions. Then for instance, multiple pions can be 
produced. If the density of the latter is sufficiently high, they can form a coherent state, thus, 
providing the feasibility of getting a pion laser [112]. When the energy of colliding nuclei 
is high, they can form fireballs of nuclear matter, having temperature and density sufficient 
for realizing the deconfinement transition [113-118]. Then, except pions, other bosons can 
be emitted, such as dibaryons and gluons, presenting the possibility of obtaining different 
nuclear-matter lasers [119]. 

In order to find out the conditions under which the massive creation of bosons in nuclear 
matter could be achieved, it is necessary to analyze the equation of state of hot and dense 
nuclear matter, taking into account various channels of reactions where composite particles 
could be formed. To our mind, the most general approach for reaching this goal is presented 
by the theory of clustering matter [117,120,121]. This approach is based on three pivotal 
concepts: cluster representation, statistical correctness, and potential scaling. 

The idea of cluster representation goes back to the works studying the abundances of 
chemical elements by treating each element as a quasiparticle characterized by its atomic 
weight and binding energy, with the related chemical potentials taking account of the allowed 
chemical reactions [122]. The problem of constructing an accurate quasiparticle representa- 
tion for composite particles was initiated by Weinberg [123-125]. The most mathematically 
elaborated approach was formulated by Girardeau [126-128], whose method was applied to 
different systems containing bound clusters, including quark-hadron matter [129,130]. 

The basic points in the cluster picture are as follows. Let us consider a multiparticle 
system with the total space of quantum states being a Fock space J-', on which the algebra of 
observables A is defined. Let among these quantum states be free states as well as different 
bound states corresponding to clusters of several particles. Each type of bound clusters can 
be individualized by a set of characteristic parameters, such as the compositeness number 
Zi, showing the number of elementary particles bound into a cluster, effective mass of the 
cluster nii, its baryon number Bi, strangeness Si, and so on. And let us treat each type of 
these bound clusters as a separate sort of particles, with the associated Fock space jFj called 
the ideal cluster space. The total cluster space is defined as the tensor product 

J^T = <^i^i®^ (201) 

termed the Fock-Tani space [126-128,131]. The relation between the Fock space of elementary 
particles and the Fock-Tani cluster space (201) is presented by means of the unitary Tani 
mapping U^ = Uj^^, so that 

jr = u^jr^ ^ jr^ = (j+jr . (202) 

The cluster algebra of observables is given by the Fock-Tani representation 

At = U^AUt . (203) 
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The latter definition guarantees that all matrix elements of the algebra Aon J-' are the same 
as those of At on J^t, since J-'j^At^t = J^^AJ-'. As far as the representation of At on J^t is 
isomorphic to that of A on J-', all observable quantities are the same in the standard picture 
of elementary particles and in the quasiparticle picture of the cluster system. 

The Tani mapping is constructed in the following way. Let aj(x) be the field operator of 
elementary particles, defined on the Fock space JF, with x being a set of spatial variables. 
Suppose ipp{xi,X2, ■ ■ ■ ,Xp) is a Schrodinger wave function describing a bound state of j9 
elementary particles. The field operator of the associated bound cluster, defined on the 
same space J-', can be presented as 

%{x) ^ J ^pix, -x,x,-x,...,xp-x) a,ix,)a,ix,) . . . «,(.,) dx.dx, ...dx,. 

The image of this bound state on the ideal cluster space jFp is given by a p-particle cluster 
described by the field operator ippix). The Tani mapping is defined by means of the unitary 
transformation 

Ut = exp (^ F 



[X 



dx 



2 

where the summation implies that over all admissible clusters. All operators ai{x) and ipp{x) 
satisfy the canonical commutation relations, depending on whether the considered particle 
or cluster is a boson or fermion. 

The algebra of observables (203) is defined on the Fock- Tani cluster space (201). In 
particular, one has the cluster Hamiltonian 

Ht = U^HUt , (204) 

where H is the initial Hamiltonian in terms of elementary particles, given on JF. The cluster 
statistical state is < At >, with the angle brackets implying the statistical averaging related 
to the cluster Hamiltonian (204). The density of the j-type cluster is 

1 - ^ f 

Pj= y < Nj > , ^J- J ^](^)^i(^) d^ ' 

where V is the system volume. We may introduce the cluster probability 

characterizing the statistical weight of the j-type clusters. This probability satisfies the 
standard properties, being semipositive and normalized, 

< Wj < 1 , y^^Wj = 1 . 
j 

The direct construction of the cluster Hamiltonian (204) is a rather tedious and com- 
plicated procedure. Moreover, the resulting Hamiltonian is presented by an infinite series. 
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Because of this it is customary to define an effective Hamiltonian whose construction is based 
on physical reasoning. Generally, an effective Hamiltonian H^/f = Heff{{pj},T) includes an 
explicit dependence on thermodynamic variables, such as cluster densities pj and tempera- 
ture T. However, including thermodynamic variables into the Hamiltonian may break the 
validity of the general thermodynamic relations. 

The correct cluster Hamiltonian has to be defined so that to preserve all thermodynamic 
relations. This is the meaning of the principle of statistical correctness [117,120]. For this 
purpose, the effective Hamiltonian is complimented by an additional term, guaranteeing the 
statistical correctness of the resulting total cluster Hamiltonian 

Ht = H,ff + CV , (206) 

which assumes the validity of the equations 

Choosing C as a nonoperator quantity yields the equations 

QC 1 ^dH,^^ dC 1 dH^ff 

defining C = C{{pj},T). These conditions guarantee the correctness of all thermodynamic 
relations, such as 

_dn __n 

~~dV~~ V ' 

dP ^ 1 rr dP I 

— -P + PbUb + Psns = - <Ht> , ^=Qf = f 



T —-- P + pbtib + Psns = — < Ht > , s = —- = —{e + P- pbtib - Psris) , 



dP dP 

nB = j— = Y. BjPj , ns = jr— = Y. SjPj ' (209) 

OPB j OPS j 

in which Q is the grand potential, e and s are the energy and entropy densities, P is pressure, 
Pb is baryon potential, ps is strangeness potential, ub and ns are the baryon and strangeness 
densities, and 

n = -T\n Tre"^^ , H = Ht-J2 A^i^i ' 

j 

dP 1 

pj = -Q— = — <Nj> , Pj = pbBj + psSj , 

with Pj being the chemical potential and [3, inverse temperature. 

The cluster Hamiltonian (206) contains the terms describing effective interactions be- 
tween different clusters. Defining the related interaction potentials is done by means of the 
principle oi potential scaling [117,120]. According to this, the interaction potentials from the 
same class of universality are connected by the scaling relation 

-?^ = ^ . (210) 

Zj^Zj ZdZf, 
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This principle allows the definition of all qualitatively equivalent interaction potentials 
through one known universal potential. Another admissible form of the potential scaling 
could be 

miirij TnaiTLiy 

since masses are usually proportional to the corresponding compositeness numbers, itlj ~ Zj. 

The theory of clustering matter has been applied to clustering quark-hadron matter, 
with elementary particles being quarks, antiquarks, and gluons, while hadrons being various 
bound clusters of these elementary particles [117,120]. 

The results of this consideration are presented in the following figures, where the energy is 
measured in units of J = 225 MeV, temperature = /c^T is given in MeV, n^B = 0.167 fm^^ 
is the normal baryon density. These results demonstrate that the deconfinement transition at 
finite baryon density, and at conditions typical of heavy-ion collisions, is a gradual crossover. 
The deconfinement transition can be associated with a point where some reduced quantities 
display a maximum. The largest amount of pions is produced around the deconfinement 
temperature ©^ = 160 MeV and at the low baryon density hb < uqb- A high concentration 
of dibaryons can be achieved at low temperature < 20 MeV and the baryon density 
ub ~ (5 — 20)noB. Dibaryons can form the Bose-Einstein condensate, thus, being in a 
coherent state. Producing high concentrations of pions or dibaryons, under the corresponding 
conditions, can be employed for creating pion and dibaryon lasers [119]. For producing a 
large amount of gluons, very high temperatures, > 160 MeV, are required. But gluons 
cannot be emitted as free particles. 

12 Conclusion 

We presented a detailed description of coherent radiation by nuclear spins. It is important 
to stress that pure spin superradiance was first experimentally observed in Dubna. The 
theoretical consideration of coherent nuclear spin radiation is based on the theory developed 
by the authors. The mathematical basis of this theory is the scale separation approach. 

We show that at sufficiently strong external magnetic fields in the system of nuclear spins 
there appear nuclear spin waves, which play the role of the triggering mechanism starting 
the incoherent motion of spins. Neither magnetodipole radiation nor the resonator Nyquist 
noise can initiate spin motion. After the spins started moving, the transverse transition 
coherence of their evolution is due to the action of the feedback magnetic field formed by the 
resonant electric circuit. This resonator field is crucial for developing coherent spin motion. 
Different regimes of coherent spin radiation are investigated. Electron-nuclear hyperfine 
coupling can lead to an essential enhancement of nuclear spin radiation. Superradiance by 
magnetic molecules is feasible, which necessarily requires the presence of a resonant electric 
circuit and cannot be achieved without it. The time dependence of the transition frequency, 
caused by the single-site magnetic anisotropy, can be compensated by the chirping effect. 
The possibility of pion and dibaryon radiation from excited nuclear matter is discussed. 

The mathematical techniques, presented in this review for treating strongly nonequilib- 
rium spin systems, are rather general and can be employed for analysing nonlinear dynamics 
of arbitrary spin or quasispin essemblies. A very interesting application of these techniques 
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would be to investigating the nonequilibrium phenomena occurring in dilute gases of cold 
trapped atoms [95,132-135]. Such atoms can be spin-polarized, forming at low temperatures 
the spinor Bose-Einstein condensates. The nonequilibrium dynamics of polarized spinor 
condensates can exhibit coherent phenomena similar to those happening in nuclear spin 

systems. 
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Figure Captions 

Fig. 1. Typical setup in experiments on spin superradiance. The left coil plays the role 
of antenna. The right coil, surrounding the studied sample, is part of the resonant electric 
circuit. The sample is shown as a dark bar inside the resonant coil. The upper left block is 
an oscilloscope, and the lower one is a plotter. The dashed line symbolizes refrigerator. 

Fig. 2. Voltage signal of a superradiant pulse as a function of time, measured in units 
of 10^^ s, for two initial spin polarizations: Sq = 0.52 (lower curve) and Sq = 0.57 (upper 
curve). 

Fig. 3. Intensity of radiation (upper curve) in arbitrary units and the longitudinal 
spin polarization (lower curve) as functions of time, measured in units of T2, obtained from 
computer simulations for 300 spins. 

Fig. 4. Orientation of the coordinate axes with respect to the spin sample inserted into 
the coil of an electric circuit. 

Fig. 5. Pressure (in units of J^) of the clustering quark-hadron matter. 

Fig. 6. Energy density (in units of J^) on the temperature-baryon density plane. 

Fig. 7. Pressure-to-energy density ratio defining the effective sound velocity squared 
clff = P/e, with a maximum at Qd = 160 MeV related to the deconfinement. 

Fig. 8. Specific heat Cy = de/dT (in units of J^). 

Fig. 9. Reduced specific heat ay = TCy/e displaying a maximum at O^ = 160 MeV 
associated with the deconfinement crossover. 

Fig. 10. Compression modulus k^^ = nsOP/dnB (in units of J^) also having a maxi- 
mum at the deconfinement crossover. 

Fig. 11. Cluster probability of the quark-gluon plasma, being the sum of the probabili- 
ties of quarks, antiquarks, and gluons. 

Fig. 12. Pion probability, being the sum of the probabilities of 7r+, tt", and n^ mesons, 
with a sharp maximum at the deconfinement crossover. 

Fig. 13. Summary probability of rj, p+, p~, p^, and uj mesons. 

Fig. 14. Nucleon probability, being the sum of the probabilities of neutrons, protons, 
antineutrons, and antiprotons. 

Fig. 15. Dibaryon probability, which is the sum of the dibaryon and antidibaryon 
weights. 

Fig. 16. Probability of the Bose-Einstein condensate of dibaryons. 
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